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ABSTRACT 


Research  on  sonic  fatigue  in  combined  environment  is  described.  Emphasis  is  placed  on 
determining  the  effects  of  structural  curvature,  low-frequency  vibratory  loads,  and  heat, 
both  singly  and  collectively,  on  sonic  fatigue.  The  analytical  and  experimental  investi¬ 
gation  is  presented  in  two  major  phases: 

1 .  An  investigation  of  simple  structural  panels  to  determine  the  effect  of 
curvature  and  heat  cycling  schemes  on  dynamic  response  and  fatigue. 

2.  An  investigation  to  determine  the  effects  of  high-intensity  sound,  heat, 
and  low-frequency  vibratory  loads  on  curved  titanium-faced  honeycomb 
sandwich  panels . 
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INTRODUCTION 


From  the  time  that  fatigue  of  flight  vehicle  structure  was  determined  to  be  partially  due  to 
power  plant  generated  acoustical  noise,  the  analytical  and  experimental  investigations 
have  been  mainly  concerned  with  flat  structure  exposed  only  to  acoustical  excitation i.  It 
Is  logical  that  this  approach  was  taken  because  sonic  fatigue  had  been  generally  confined 
to  secondary  structure,  control  surfaces,  and  fairings.  These  structural  components  are 
inherently  flat  and  generally  exposed  to  the  most  severe  acoustical  noise  *  j*  * 

during  on-ground  operation  and  high  gross  weight  takeoff.  During  these  segments  of  flight 
vehicle  operation  the  air-loads  and  aerodynamic  effects  are  minimal. 

The  assumption  that  the  useful  life  of  a  structural  component  can  be  determined  from 
consideration  of  acoustical  excitation  alone  is  one  that  is  open  to  considerable  doubt .  For 
example,  during  the  sonic  proof  test  of  the  B-58  Weapon  System  it  was  determined  con¬ 
clusive!  v  that  low-frequency,  wing-mode  vibratory  loads  caused  a  large  percentage  ot  the 
fatigue  damage  in  certain  wing  structure  (Reference  1).  Body-mode  stresses  are  certain  y 
more  intense  during  taxi  and  high  speed  takeoff.  It  is  quite  possible  that  these  dyrwmic 
loads  have  caused  some  of  the  so  called  "early"  sonic  fatigue  failures  of  flight  vehicle 
structure  subjected  to  both  types  of  loads.  In  addition  to  the  multiple  loads  creating 
fatigue,  some  of  the  structure  is  immersed  in  the.high  thermal  environment  of  the  let  eff  ux 
with  an  attendant  reduction  in  the  allowable  fatigue  strength. 


High  performance  flight  vehicles  will  experience  fatigue,  significant  boundary  layer 
turbulence,  aerodynamic  heating,  and  gust  or  buffet  effects  simultaneous  y .  ®SIPI1 
criteria  for  advanced  flight  vehicles  are  to  be  realistic,  these  environmental  conditions 
can  no  longer  be  separated  as  they  have  in  the  past* 

To  this  end,  the  research  reported  herein  was  designed  to  determine  the  effects  of  structure 
curvature,  low-frequency  vibratory  loads,  and  heat,  both  singly  and  collectively,  on  some 
fatigue.  The  overall  investigation  consisted  of  two  major  phases.  These  were: 

1  .  An  investigation  of  simple  structural  panels  to  determine  the  effect  of  curvature 
and  heat  cycling  schemes  on  dynamic  response  and  fatigue. 

2.  An  Investigation  to  determine  the  effects  of  high  intensity  sound,  heat,  and 

low-frequency  vibratory  loads  on  curved,  complex  structural  panels  representative 
of  future  flight  vehicle  structure. 


A  description  of  the  Investigations  and  results  of  the  first  phase  is  contained  in  Section 
The  investigations  and  results  for  the  complex  structural  panels,  the  second  phase,  are 
described  in  Section  III. 


1 


II 

SIMPLE  SPECIMEN  INVESTIGATION 


A.  Theoretical  and  Analytical 


1 .  Statement  of  Problem 

The  method  chosen  for  defining  the  basic  effects  of  curvature  on  fatigue  life  is  to  test  to 
failure  the  most  simple  curved  panel  specimen  that  could  be  designed.  This  configuration 
is  an  unstiffened,  cylindrical  ly  curved  panel  which  is  fastened  to  a  semi-rigid  frame. 

This  section  of  the  report  is  a  description  of  the  analysis  which  defines  the  panel  response 
to  an  arbitrary  acoustic  excitation. 

Three  of  the  basic  approaches  which  can  be  pursued  in  a  response  analysis  of  an  elastic 
structure  are: 

a.  the  differential  equations  approach 

b.  the  lumped-mass  analogy,  and 

c.  the  energy  approach 

The  differential  equations  approach  is  generally  termed  an  "exact  analysis"  within  the 
limits  of  linear  homogeneous  elasticity  relationships  and  other  simplifying  assumptions 
which  are  used.  The  differential  equations  for  a  complex  structural  system  can  be 
developed  by  equilibrium  considerations.  Also,  a  variational  approach,  utilizing  the 
Euler-Lagrange  equations  in  conjunction  with  the  total  energy,  can  be  used  to  derive  the 
differential  equations.  It  is  exceedingly  difficult,  however,  to  obtain  solutions  to  the 
differential  equations  for  other  than  simply-supported  boundary  atmeHLticafj  unless  approxi¬ 
mate  methods  such  as  a  finite  difference  or  Gallerkin  are  utilized. 

The  lumped-mass  analogy  is  also  straightforward  to  set  up  and  solve.  This  method  can 
handle  most  boundary  conditions,  but  it  is  useful  only  If  a  digital  computer  is  available. 

It  does  not  rely  upon  assumed  mode  shapes,  but  uses  standard  methods  for  generating  the 
lumped-mass  system  modes.  The  one  disadvantage  of  this  method  is  that  when  a  two- 
dimensional,  or  lattice-work  structure  is  analyzed,  even  the  largest  computers  do  not 
have  adequate  capacity  to  solve  the  problem,  considering  the  number  of  mass  stations 
necessary  to  adequately  represent  the  system. 

The  Rayleigh-Ritz  energy  method  is  exact  within  the  limitations  of  the  differential 
equations  approach  with  one  exception:  a  mode  shape  must  be  assumed  which  is  not 
necessarily  a  solution  to  the  differential  equations.  In  cases  where  exact  mode  shapes  can 
be  assumed,  the  final  results  are  exact.  For  boundaries  such  as  clamped  edges,  inexact 
mode  functions  which  satisfy  only  the  geometric  boundary  conditions  can  be  used.  The 
Rayleigh-Ritz  energy  solutions  have  been  relied  upon,  mainly  because  results  are 
guaranteed . 

The  basic  approach  to  the  curved  panel  response  has  been  that  of  solving  for  both  the 
simply-supported  and  clamped  edge  conditions.  Since  most  structural  panels  have  edge 
fixities  which  lie  between  these  boundary  conditions,  the  two  analyses  provide  an  upper 
and  lower  bound  for  both  frequency  and  response. 
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Since  the  modal-energy  approach  has  been  used,  the  problems  of  defining  panel  response 
logically  divides  into  two  separate  problems,  i.e.  (a)  definition  of  natural  frequencies, 
mode  shapes,  and  generalized  mass,  and  (b)  definition  of  response.  These  topics  comprise 
the  following  two  sections  of  the  report. 


2,  Frequency  Analysis 

The  vibration  analysis  of  a  cylindrical ly  curved  panel  was  made  using  the  following 
assumptions: 

a.  The  material  is  linearly  elastic. 

b.  The  material  is  orthotropic  . 

c.  Panel  thickness  is  much  less  than  the  major  panel  dimensions.  (Elasticity  theory 
of  thin  shells  is  applicable.) 

For  ease  in  accomplishing  the  manipulations,  matrix  algebra  will  be  used  throughout  the 
analysis . 

The  strain  energy  density  of  an  elastic  system  is 

UQ  =  2  LPjJ  C e 1 1 


It  is  necessary  to  express  the  stress,  a*,  in  terms  of  strain,  e.,  and  then  the  strain  in  terms 
of  displacements .  For  an  orthotropic  1  elastic  solid,  the  stress-strain  relationships 


are 


{0.3  *  [c..]  Uj)  LCT[J  ~  LcjJ  C.J  (2) 

where  IC..J  is  the  elastic  coefficient  matrix.  For  orthotropic  materials  C  is  symmetric. 

1 1 

The  strain  displacement  relationships,  as  derived  by  Love  (Reference  2)  ore 
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ay 

a 

< 

>  -  z 

a 

ay2 

< 

v> 

a 

a 

0 

w 

0  lit 

a2 

axay 

— 

*  J 

This  equation  can  be  written  in  shorthand  form  as 

[e.)  =  [[p]  -z  [Q|]  <v  ► 

[wj 

u,  v,  and  w  are  shown  in  Figure  1,  which  defines  the  coordinate  system. 


(3) 


(4) 
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FIGURE  1.  CURVED  PANEL  COORDINATE  SYSTEM 
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The  boundary  conditions  for  a  curved  plate  with  clamped  edges  are 

w(0,y)  =  w(£,y)  =  w{x,  0)  =  w(x,b)  =0 

wx(0,y)  =  w^,y)  =  w  (x,0)  =  w  <x,b)  =  0 

v(0,y)  =  v(t,y)  =v(x,0)  =  vjc,b)  =0 

u(0,y)  =  u(t,y)  =  u(x,0)  =  u(x,b)  =0 

The  geometric  boundary  conditions  for  a  curved  plate  with  simply-supported  edges  are 

w(0,y)  =w(t,y)  =  w(x,0)  =  w(x,b)  =  0 
vy(x,0)=vy(x,b)  =  0 
ux<°,y)  =  ux(-t,  y)  =  0 

Beam  deflection  functions  of  the  form  shown  below  satisfy  both  boundary  conditions. 


u  'II 1“  Umn  XmW  YnW  1 

X,M 

*  "ll  W„n  Xn,«  YnW 


For  clamped  edges 


X  (x)  -  Cosh  0  x  -  Cos  0  x-a.  (Sinh  0  x  -  Sin  0  x) 
m  m  mm  m  m 

Yn(y)  =  Cosh  Yny  -  Cos  Yny  "  9n(Sinh  Yny  "  Sin  Yny) 


a  = 


Cosh  p  t  -  Cos  p  £ 
_ _ m _ m 

m  Sinh  p  t  -  Sin  p  t 
m  m 

Cosh  y  b  -  Cos  y  b 
a  =  n _ ln_ 

n  Sinh  y  b  -  Sin  y  b 


and  and  Yn  are  determined  from  the  following  relationships: 


Cosh  P  t  Cos  P  l  =  1 
m  m 

Cosh  y  b  Cos  y  b  =  1 
n  Tn 


The  assumed  deflection  functions  for  simply -supported  edges  are 

X  (x)  =  Sin  P  x 


m 


m 


Yn(y)  =SinYny 


(5) 


(6) 


(7) 


(8) 


Equation  (5)  in  matrix  form  is 
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(9) 


0  0 


Substitution  of  (9)  into  (4)  yields  the  strain-d 
shapes  and  generalized  coordinates, 


- 

p  ■> 

0 

Umn 

.,1 

/mn 

0 

< 

V 

mn 

it 

/> 

V  ^ 
mn 

LXm  V. 

W 

mn 

I*  * 

W 

mn 

t-  J 

placement  relationships  in  terms  of  the  mode 


[e,l  -[|>]-*[Q]]  M 


=  [[b]-z[d]]  tqr] 


(10) 


The  transpose  of  {e,}  is 

LEjJ  =  LqsJ  [[b]T  -  z[dJT]  (11) 

Substituting  (2),  (10),  and  (11)  into  (1)  results  in 

h,4lv[[b]t-M][c][[bMd]]  h!  os 

Performing  the  above-indicated  matrix  multiplications  yields 

^lv[[b]tMM-M[c][b] 

-,[bJMc][d]  +  ,2[d]t  [c][d]  ]ur]  (13) 

The  total  strain  energy  of  the  curved  plate  is  obtained  by  integrating  the  strain  energy 
density,  U  ,  over  the  volume  of  the  plate.  In  symbolic  form,  the  strain  energy  is 

b  t  h/2 

U  ~Jff  UQ  dz  dx  dy  (14) 

0  6  -h/2 

Integration  across  the  thickness  gives 

u=f  f  /  LqsJ  /Mt[c][b]  +  ^[d]t[c:][d]]  (qr)  dxdy  05) 

6  o 

Following  through  with  the  indicated  matrix  multiplications  results  in 
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(16) 


U  = 


iff  LqsJ  [h[E]+^[F]]{qr}dxdy 
0  0 


(Note:  The  terms  contained  in  [e]  are  the  stretching  or  membrane  contribution  to  strain 
energy  whereas  [Fj  constitutes  the  bending  energy.)  Integrating  over  the  surface  area  of 
the  plate  yields 


(17) 


The  kinetic  energy  of  the  vibrating  plate  is  obtained  by  Integrating  the  product  of  mass 
and  one-half  velocity  squared  over  the  volume  of  the  plate: 

b  6  h/2 

T  =  £  J  j  J  (u^  +  v^+\v^)dzdydx  (18) 

0  0  -h/2 

Substituting  (9)  for  displacements  into  (18)  and  performing  the  indicated  time  derivatives 
produces 


b  l  h/2 

H///  L6mnVmnWmnJ[AjT[A] 

0  0  -h/2 


U. 


mn 


,mn 


W 


mn 


i  dz  dx  dy 


(19) 


Integration  over  the  volume  results  in 


T=^[j](qr] 


Lagrange’s  equation  of  motion  is  used  to  obtain  natural  modes  of  vibration  as  follows: 

I1  -  li<J  (qr) 


£(J£)=PHUJ  fq'r} 


(20) 

(21) 

(22) 


Assuming  harmonic  motion,  (22)  becomes 

lr(^)=-uj2phLJJ  Hi 


(23) 


Substituting  (21)  and  (23)  into  LaGrange's  equation  results  in  an  eigenvector  solution: 

[  [<]  -  uu2  ph[j  ]]  Eqp  =  0  (24) 
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The  intermediate  manipulations  required  to  derive  the  terms  in  the^X]  and[j]  matrix  will 
not  be  presented  for  the  sake  of  brevity.  The  terms  In  the  fx]  matrix  can  be  partitioned 
into  a  3  x  3  symmetrical  matrix.  The  elements  of  this  matrix  are  square  matrices  of  order 
MN,  These  terms  are 


Kn=h{ 

K12=  h  {C 
K 1 3  =  h  {c 


cn  A 


3M 


pm 


3  V 

pm 


JNqnJ 


'66 


2m 


pm 


P  V 

urp  m 


[2Nqn]  } 


12 


> 


pm 


L 


3  t 


[§3]*^  [*=&]]} 


',M  r 

'12  ^L"^1  JNqn 


■]]> 


q  n 


+  C66A(,+^) 


ijM  I 

l  pm 


&]]> 


>(25) 


y  Y  b 

q  " 


K23 


Sr  S. 

KA 

2Nqn 

.Si* 

KA 

i 

hO 

Z 

_Q 

13 

_ I 

" 

1^22  A 

1  tv\ 

L 1  Pm 

LvqM 

.  12At2 

ajVI 

l  pm 

y  b 

L  Tq  J 

C22<P 

1 2t2A3 


M*]] 


c66t<T 


3A7L  A" 


[*]1) 


K33  h{C22  AT  [lMpm  [lNqn]J  +  [3Mpm  [  1  Nqn]] 


'22 


1 2t2A3 


[lMpm[  3^qn]j  +  <2C12  +  4C66*  ]  ^  [2Mpm[2Nqn]_  } 


The  [  j]  matrix  can  be  shown  to  be 

J, ,  = 


J22 


[l^T[b  lNqJ 

=  lb  | 

pm 

1 

] 

[* >M4$£] 

-  lb  j 

j 

[  '“-01 
q  n 

r  t 

]  ’ 

Ct  ,M  f  b  ,N 

1  pmL  1  qn, 

]  =  tb 

LlMpm[lNqn]J 

(26) 


The  non-dimensional  parameters  A  =  b/l,  cp  =  b/a,  and  t=-t/h  were  used  inequations  (25)  and  (26) 
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Since  the  matrix  notion  j^Mpm  [Nqn]  j  is  unconventional  the  following  definition  is  given: 


M  = 


mll  m12 
m21  m22 


[n]  - 


M  - 


m^Ll]  mi  2^l] 

m2-|[l] 

[  1  ]  is  the  identity  matrix  of  the  same  order  as  [N] 


M  ■ 


nn  ni2 


n21  n22 
[Nj  0 
0  [Nj 


[m[n]]  =  [m]  [n1]  = 


m^  [N]  m12  [N] 

[Nj  m22  f  N] 


L  21 


Integrations  of  the  clamped  mode  functions  over  the  surface  which  were  used  to  derive  (25) 
and  (26)  are  (Reference  3) 


X  X  dx  =  6  t 

J  p  m  pm 


f  X1  X'  dx  =  6  a  pja  p  t  -  2) 

J  pm  pm  m  m  m  m 


=  t  ,M 

1  pm 


+  o-6  )  4P^(y m ;  w  [,  +  H)^p],  (v 


pm 


4  4 

p  -  P 

m  p 


f  X"  X"  dx  =  6  B*  t 
J  pm  pm  m 

o 

i  t 

f  X"  X  dx  =  f  X  X"  dx 
J  p  m  J  P  m 


f  V  Y  dy  =  6  b 
J  an7  qn 


-M 

_  3  pm 

" 


L  pm 

l 


=  b  ,N 

1  qn 


(27) 
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f  Y1  Y1  dy  =  6  0  y  (9  y  b  -  2) 

J  q  n  7  qn  n  T  rV  nTn 


4y^  Y^(®  Y  -  0  Y  )  r 

+  0  -  6J  q  Vq4  ""  ['  +  H)^"]  =  V2 

H  Y„  -  Y, 


'n  q 


I  Y"  Y"  dy  =  6  y4  b 
J  q  n  7  qn  n 


0 


b  b 

f  Y  YM  dy  =  f  Y"  Y  dy 
J  q  n  7  J  q  n  7 

0  0 


where  6..  Is  the  Kroneker  delta: 
<1 


^  (27) 

3N  (cart'd) 


«N 

_  l  qn 

b 


6..  = 

'I 


1  i  =  j 

0  i^j 


Table  1  is  a  tabulation  of  the  non-dimensional  values  ..  -  2^j"  These  values  are  taken 
from  Young  (Reference  4). 


TABLE  1 

VALUES  OF  ~M..  =  „N.. 

2  i I  2  i) 


i/i 

1 

2 

3 

4 

5 

6 

1 

12.30262 

0 

-9.73079 

0 

-7.61544 

0 

2 

0 

46.05012 

0 

-17.12892 

0 

-15.19457 

3 

-9.73079 

0 

98.90480 

0 

-24.34987 

0 

4 

0 

-17.12892 

0 

171  .58566 

0 

-31 .27645 

5 

-7.61544 

0 

-24.34987 

0 

263.99798 

0 

6 

0 

-15.19457 

0 

-31.27645 

0 

376.15008 

Table  2  is  a  listing  of  the  non-dimensional  quantities  PX  -  y-b  taken  from  results  published 
by  Young  and  Felgar  (Reference  5). 
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TABLE  2 

VALUES  OF  £X  “  y.b  FOR  CLAMPED  BOUNDARIES 


i 

BX  =  y.b 

(B.t)4  -  (y.b)4 

1 

4.7300408 

500.5639 

2 

7.85320046 

3,800.5370 

3 

10.9956078 

14,617.6299 

4 

14.1371655 

39,943.7991 

5 

17.2787596 

89,135.4065 

6 

20.4203522 

173,881 .316 

Clamped  Boundaries:  i  >  6;  0X  -  y.b  -  (2i  +  1)  ^ 

Simply  Supported  Boundaries:  i  =  1 , 2, . . 

8X  =  y.b  -  in 


The  values  of  the  integrals  of  the  assumed  mode  functions  over  the  plate  surface  for 
simply -supported  edge  conditions  are 

l 

f  X  X  dx 
J  pm 


0 


f  X"  X  dx  = 
J  pm 


f  Y  Y  dy  =  6  fe 
J  q  n  /  qn  / 


t 

)  rt 

pm  l 

=  6  ,M 

1  pm 

s  g2  4 

pm  ml 

0M 

2  pm 

-t 

> 

pm  m  l 

-M 

■V 

l 

f  X  X" 

J  pm 

D 

-M 

b 

qn  2 

=  b  ,N 

1  qn 

2  u 

Y  b 
n 

*qn 

,N 

l  qn 

b 

(28) 
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b  4,  M 

c  v  b  ,N 

f  Y"  Y"  dy  =  6  -4-  =  t11 

J  q  n  '  qn  l  Lo 

0  b 

b  b 


9N 
l  cjn 


0  0 

j 

For  an  orthotropic  elastic,  thin  plate,  the  elastic  coefficient  matrix  reduces  to 

0 


W- 


where 


C11  C12 

C21  C22 
0  0 

E 

C  = - - — 

^11  1  -u  u 

x  y 

u  E 

r  = _ y  y 

1 2  1  -u  u 

x  y 

u  E 

r  = _ v  x 

*“21  1  -  u  u 

x  y 


'66 


C„  =  t - * - 

22  1  -  u  u 

x  y 

C,,  =  G 
66  xy 

and  the  relationship  ux  Ey  =  Uy  Ex  is  necessary  for  orthotropic  media,  and  the  elastic 
coefficient  matrix  must  be  symmetric  (Reference  7). 


(28) 

(cont'd) 


(29) 


This  concludes  the  frequency  analysis  for  orthotropic  curved  panels;  however,  there  are 
several  important  points  which  need  to  be  discussed.  These  are 

a.  Referring  to  equation  (27),  the  integrals  of  Xp  X^,,  Xp  Xm  and  XpX^,  for 
clamped  edge  conditions  are  non-zero  when  p  is  not  equal  to  m.  Thus  the 
analysis  doesn't  display  the  desired  orthogonality  between  the  modes.  A 
numerical  analysis  for  one  of  the  test  panels  used  in  this  program  compared 
with  a  numerical  analysis  which  assumed  orthogonality  shows  insignificant 
differences;  however,  a  complete  investigation  of  the  effects  of  including 
this  non-orthogonality  relationship  has  not  been  evaluated  because  of  com¬ 
puter  time  requirements.  In  order  to  obtain  the  first  25  (through  5,5  mode) 
panel  modes  from  (24)  for  clamped  edges,  it  is  necessary  to  solve  a  75  x  75 
eigenvalue  matrix.  Computer  time  required  for  this  analysis  precluded 
parameter  studies . 
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b.  Examination  of  (28)  reveals  that  the  integrals  of  the  mode  shapes  are 
orthogonal  for  simply-supported  edges,  since  sine  functions  were  chosen. 

c.  A  simplification  of  considerable  interest  to  the  orthotropic  curved  panel 
frequency  analysis  occurs  provided  the  modal  integrations  are  taken  to  be 
orthogonal  and  the  material  is  isotropic.  In  this  case  [iMpnJ,  [^MpmJf 
[3Mpm],  ClNqn],  [2Nan]  and  [jjNqn]  become  diagonal  matrices  and  the 
modes  are  uncoupled.  Also,  the  elastic  matrix  reduces  to 


C11  ”  C22 


1  -u' 


r  -  r  —  V  E  _ 

l12  ^21  ,  2 

1  -  u 


C66 =  G  = 


2(1  +  v) 


Also,  it  is  assumed  that 


:«  1 


Incorporation  of  these  parameters  into  the  eigenvalue  problem  results  in  a  3  x 
3  matrix  for  determining  the  m  n*b  natural  frequency.  The  determinant  of 
coefficients  is 


W-x2W|-o 


where 


H 


and 


L  i*(l  -  u  ) 

[  j]  =  -tb  [l] 


^2=p<.3  b(l  -u2)^ 
Efi 


The  terms  in  the  G  matrix  are 

-  _  .2  f]  +  u\  M2(m)N2(n) 

G12~'  CT7  (B  t)(v  b) 


G11  =  At 


2M3(m)N1(n)  f2  ,  _  .  M2(m)N2(n) 

I'/  +  A  ^  2 


m  '  n 


(30) 
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2  M2(m)Nj(n) 

G 1  j  =  cpt  vj 


p  t 
rm 


i2  «  /i  _  N2(n) 

G22  =  A  M1M  .  ' ,  v2  +  At  (  2  )  M2^  ,,,5 
(y  b)  (Y  b) 


.2  N»  /0  »  N» 

G23=^M1^^T+|(^M2H— + 


T  ’  .  <P  tJ  ,  ,  3 


Y  b 


~  M,  (m) - r 

12A3  1  Ynb 


N,(n) 


=  —  *  (m)N^  (n)  +  ^ 2  M3(m)N  j  (n)  +  — ^  M  j  (mjN^n)  +  A/^mjNjtn) 

A  1 2A  6A 

The  terms  in  the  matrix  are 


LI)  = 


M2(m)  N|{n) 

0  i)2 

m  ' 

N,(n) 

L22  ~  M|(m>  2 

un  / 

L33  =  M](m)  N^n) 


This  form  of  the  frequency  theory  reduces  to  an  analysis  performed  by  Sewell 
(Reference  6) . 

Table  3  contains  a  tabulation  of  M;  and  N;  for  both  clamped  and 
simply-supported  edges. 


(30) 

(cont'd) 


TABLE  3 

VALUES  OF  M.  AND  N. 

i  i 


Clamped 

Simply-Supported 

Mj  (m) 

i 

1/2 

M2(m) 

am  ^  ^  St  "  2) 

Hi  m  mm 

(3m£)2/2 

M3(m) 

(e  D4 

m 

®J)4/2 

N,(n) 

1 

1/2 

N2(n) 

0  y  b(9  y  b  -  2) 

(Y„b)2/2 

N3(n) 

(Ynb)4 

(Ynb)4/2 
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d .  Finally,  if  a  =  »  (flat  plate,  cp  =  0)  then  the  3  x  3  matrix  reduces  to  a  2  x  2 
and  one  equation  in  terms  of  X2  in  the  3,3  position.  The  equation  resulting 
from  the  3,3  element  yields  the  flat  plate  flexural  modes,  whereas  the  2  x  2 
matrix  gives  the  in~plane  or  longitudinal  vibration  modes. 


3.  Response  to  Acoustical  Excitations 

The  analytical  procedure  used  for  evaluating  panel  response  to  an  arbitrary  acoustic 
excitation  is  a  modal  analysis  derived  from  LaGrange's  equations  having  the  form 
(Reference  8): 


d  pj.'i.jji.+aJi +JS..Q. 

dt  \dq./  dq.  dq.  dq.  ■ 


s.Sq  j  /  dq.  dq.  dq. 

where  R  is  the  dissipation  function  defined  as 

R  =  2pT  ,  8  =  (c/c  )(ju. 

C  I 

and  it  is  necessary  to  determine  dT/dq.,  dT/dq;,  and  dU/dq.. 


(31) 


(32) 


Although  the  modes  do  not  display  orthogonality  between  some  of  the  terms  (see  (27)), 
orthogonality  will  be  assumed  for  the  response  calculation.  From  (22), 


dq- 


phL-Jj  -  tqj} 


(33) 


However,  since 


[j]  is  diagonal, 


also 


#=phj.q.  =  M.  q- 
dq.  i  ~i  i  ' i 


f-  =  o 

dq. 


(34) 


|^-=  LKJ,  tq:)  (35) 

dq;  >  i 

If  the  system  is  vibrating  in  the  i ^  natural  mode,  with  zero  damping  and  harmonic  motion, 
(31)  becomes 


dU 

dq. 


+ 


±{$L)  = 

dt  \dq./ 


[LKJ; 


phLJj.]  {qj} 


(36) 


For  arbitrary  values  of  q.. 


2  2  . 
LKJ  .  =  uu.  ph|_Jj  .  =  uu.  M* 


(37) 
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Substituting  (32),  (34),  and  (37)  into  (31)  produces 


M.  q.  +  — a).  M.  i\.  +  uu.  M.  q.  =  Q. 
i  ~i  c  i  i  ~l  1  I  I 

c 


(38) 


For  a  distributed  acoustical  pressure,  Q.  has  the  following  form: 

b  t 

Qj  =  /  /  p(x/y/t)(X(x)Y(y)).  dxdy 


(39) 


0  0 


p(>t/y/0  is  the  acoustic  pressure  distribution. 


th 


For  the  mn  mode,  using  the  results  of  (30),  (38)  reduces  to 

b  l 

M  W  +  uu  M  W  +  iD^  MW  s  f  I  p(x,y,t)X  (x)Y  (y)dxdy  (40) 
mn  mn  c  mn  mn  mn  mn  mn  mn  j  J  '7  '  '  mv  7  n  7  7  '  ' 

c  0  0 

where  Mmn  and  (Dmn  are  the  generalized  mass  and  natural  frequency  associated  with 
resonance  in  the  radial  direction  (i.e.,  transverse  vibration). 


The  solution  of  (40),  assuming  harmonic  motion,  is 

Q. 


W  = 


mn 


mn  M  a.2 

mn  mn 


- 

/ 

,2 

/  \2l 

■ 

1  -  JU—\  +  , 

( or  at  \ 

\ 

'  at  )  +' 

^  ^mn  id  / 

mn 

mn  J 

1-1/2 


S  i  n(uut  +  cp  ) 


i  r  c  (uu/uu  ) 

-t  smn  '  '  mn 


CD  =  tan  « 

i  -  m°.  )2 


mn 


b  l 


Qmn  "/  /  P<x^'t)Xm(x)Yn(y)dxdy 


0  0 


For  uniform  pressure  distribution,  Qmn  reduces  to 


a  0 


(41) 


Qmn  =  4f‘rr![,-(-,)m][l-H)n] 
m  n 


for  clamped  edges  and 


m  n 


for  simply  supported  edges. 


(42) 
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(43) 


The  generalized  mass  associated  with  radial  vibration  is 

Mmn  =  pht-b  for  clamped  edges 

Mfnn  =  ^4^  ^or  simply  suPP°rtecl  edges 

Upon  definition  of  Wmn  and  uumn,  Umn  and  Vmn  can  be  determined  from  (30). 

A  procedure  is  now  available  for  determinating  the  absolute  magnitude  of  Umn,  Vmn  and 
Wmn,  The  displacements  u,  v,  and  w,  can  be  determined  using  (5). 

The  surface  stresses  for  a  thin  cylindrical  shell  are  (Reference  9) 


In  the  evaluation  of  the  effects  of  curvature  upon  natural  frequencies,  stress  response  and 
ultimately  fatigue  life,  it  is  desirable  to  arrange  the  theory  to  yield  frequency  ratio  and 
stress  ratio . 

Some  important  simplifications  can  be  made  in  the  frequency  theory  if  the  angle  which  the 
panel  subtends  is  small.  For  angles,  cp,  less  than  0.2  radians,  the  frequency  of  flexural 
vibration  can  be  approximated  by 


X  K33//J33 


(45) 


With  all  edges  clamped,  this  reduces  to 

.2  .  25.2  ,  41  .7  ,  *V 

X  =  41 .7A  +  ^  •  +  — 2"  + - 

A  A 

Equations  (44)  for  maximum  stress  on  the  straight  edge  reduces  to 

2 


a  = 


rav  h  awl 

L^y '  2  rrJ 


y  i  -u-  L“r  -  ay- 
Substituting  the  mode  shape  at  the  center  of  the  straight  edge  gives 

1  YnXm<^Wmn{7T  ’  '} 


a  =  — - - 2' 


y  i  - 


(4b) 


(47) 


(48) 


v 


17 


K  is  the  ratio  of  the  generalized  coordinates  Vmrv/Wrrin  and  is  expressed  by  expanding  the 
first  two  equations  which  make  up  the  frequency  determinant,  (30). 


[k2i  K13  -  K23^K)1  ~ Nrin  Jjp]  Wmn  _ 

<Kll-XmnJn)(K22-XmnJ22>-K21K12 


(49) 


For  thin  panels. 


X2J 


n 


«  K  i  ^  and 


\  J  22  ^ 


‘22- 


Using  this  approximation 


9nh 


A_ +  0.108 

A4  +  9.62A2  +  T 


Atcp 


(50) 


The  theoretical  first  mode  value  of  D  is  0.445  for  a  panel  with  clamped  edges. 


The  ratio  of  curved  panel  stress  to  flat  panel  stress  at  resonance  reduces  to,  using  results  of 
(48)  and  (41), 


2 


From  (46), 


1  + 


CfAtcnV 


A4+  0.61A2  +  1 


(52) 


The  theoretical  value  of  C  is  0.024  for  clamped  edges.  Combining  the  results  of  (50), 
(51)  and  (52)  gives 


C  (Atcp) 2 

I  +  — j — - - - 

A4+  0.6]A*+  iJ 


+  0.108 

9.62A2  +  1 


(53) 


A  design  nomograph  which  solves  (53)  is  included  in  Appendix  I  •  This  nomogram  is  to  be 
used  in  conjunction  with  existing  design  charts  for  flat  rib-skin  structure*  The  constants  C 
and  D  were  determined  from  the  experimental  test  results  discussed  in  Section  11. B* 


B.  Experimental 


I  ,  Introduction 

This  experimental  investigation  was  designed  to  determine  the  effect  of  curvature  on 
structural  response  and  acoustical  fatigue  and  the  effect  of  thermal  cycling  on  acoustical 
fatigue*  Broad-band  high  intensity  sound  tests  on  simple  panels  were  made  to  evaluate 
these  effects. 
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To  determine  the  effect  of  curvature,  simple  unstiffened  panels  were  used.  The  panel 
thickness  and  size  were  held  constant  within  manufacturing  tolerances  and  the  radius  of 
curvature  was  varied  to  encompass  a  forseeable  range  for  present  and  future  flight  vehicles. 
Flat  panels  made  of  heat  resistant  materials  were  used  to  evaluate  the  thermal  cycling 
effects . 

A  description  of  the  test  facility,  data  collection  and  data  reduction  system,  test 
procedures,  and  test  results  is  contained  in  subsequent  sections. 

2.  Laboratory  Test  Facilities 

The  High  Intensity  Sound  System,  designed  by  Ling-Altec,  Inc.,  is  an  electronically 
controlled  system  which  can  produce  random,  sinusoidal,  or  a  taped  output  as  desired. 

The  system  has  the  capability  of  producing  an  overall  sound  pressure  level  output  of  163 
decibels  when  operating  in  the  random  mode.  The  useful  frequency  range  in  the  random 
mode  is  50  to  2000  cps.  Random  mode  operation  is  adjustable  in  octave  band  widths. 

Air  for  the  system  is  supplied  from  plant  compressors  and  is  used  up  to  50  psig.  The  air 
flow  is  modulated  by  an  electrical  signal  to  the  four  transducers  in  the  system  to  produce 
the  desired  acoustical  noise  output.  The  sound  passes  through  exponential  horns  to  the 
progressive  wave  test  section  where  it  impinges  on  test  specimens  with  grazing  incidence. 
The  test  section  is  one-foot  wide  and  can  accommodate  test  specimens  4  feet  by  10  feet  in 
size.  Figure  2  is  a  photograph  which  shows  the  general  arrangement  of  the  system. 

The  complete  complement  of  sound  generating  system  controls,  data  collection,  and  data 
reduction  system  is  shown  schematically  in  Figure  3.  Specifically,  transducer  (micro- 
hones,  accelerometers,  and  strain  gage)  outputs  are  displayed  on  meters  and/or  oscil- 
oscopes  and  recorded  on  a  fourteen-track  Ampex  FM  magnetic  tape  recorder.  Data 
processing  equipment  includes  an  Ampex  tape  loop  system  with  data  analyses  performed 
using  a  Technical  Products  Corporation  Wave  Analyzer  system  with  a  filter  having  a 
nominal  bandwidth  of  one  cps.  A  Bruel  and  Kjaer  one-third  octave  band  analyzer  and 
plotter  is  used  mostly  to  reduce  acoustical  noise  data. 

3.  Ambient  Temperature  Tests 

The  simple  test  specimen  response  and  fatigue  tests  were  conducted  at  ambient  temperature 
to  minimize  the  variables.  A  description  of  the  test  specimens,  test  procedure,  and  test 
results  is  given  below: 

a.  Test  Specimen.  For  reasons  related  to  stress,  linearity,  and  time  to  failure,  two 
panel  co  nt  ig  u rat io ns  were  tested.  The  first  configuration  had  flat  pattern  overall 
dimensions  of  1 1"  x  13"  x  0.051"  with  arbitrarily  chosen  radii  of  curvature  of 
48",  72",  96",  and  infinity.  The  second  configuration  was  the  same  as  the  first 
except  for  on  0.032"  thickness.  Both  configurations  were  7075-T6  clad  aluminum 
alloy. 

Only  a  few  of  the  0.051"  test  specimens  were  tested  because  the  time  required 
to  fatigue  the  panels  with  the  48"  radius  was  unreasonable.  Consequently,  the 
0.051"  panels  were  discarded  in  lieu  of  the  0.032"  panels. 

Figure  4  is  a  drawing  which  gives  the  details  of  the  0.032"  test  specimens. 

b.  Test  Procedure .  In  conducting  the  fatigue  tests  the  following  procedure  was 
employed . 
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FIGURE  2.  GENERAL  VIEW  OF  HIGH  INTENSITY  SOUND  SYSTEM 
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FIGURE  3.  BLOCK  DIAGRAM  OF  HIGH  INTENSITY  SOUND  SYSTEM  CONTROLS,  DATA  COLLECTION,  AND  DATA 
AND  DATA  REDUCTION  SYSTEM 
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First,  mode  shapes  for  each  panel  configuration  were  determined  *  The  panels 
were  fastened  in  a  picture  frame  mounting  fixture  of  the  type  shown  in  Figure  5. 
They  were  mounted  one  at  a  time  over  a  loudspeaker.  The  modes  were  excited 
and  the  Cbladnr  patterns  were  photographed,  when  possible,  through  the  (5,5) 
mode.  After  a  study  of  the  mode  shapes  and  panel  response  characteristics,  it 
was  decided  that  strain  gages  should  be  placed  in  the  four  locations  shown  in 
Figure  6  on  one  panel  of  each  configuration.  Frequency  sweeps  were  then  made 
at  various  discrete  Frequency  sound  pressure  levels  to  determine  the  characteristics 
of  each  panel  response,  damping  ratios,  and  the  range  of  linear  response. 

Strain  gaged  specimens,  mounted  as  shown  in  Figure  7,  were  next  excited  by 
low-level  broad-band  acoustical  noise  having  a  flat  spectrum  to  determine  the 
most  active  strain  gages  and  the  significant  modes  of  response.  The  data  obtained 
from  the  sinusoidal  frequency  sweeps  and  the  low-level  broad-band  excitation 
were  used  to  establish  the  most  active  strain  gages  to  be  used  for  monitoring  pur¬ 
poses  during  the  fatigue  tests  and  to  shape  the  broad -band  acoustical  noise 
spectra . 

Last,  fatigue  tests  were  conducted  using  the  broad-band  spectra  shaped  from 
frequency  sweep  and  low  level  broad-band  excitation  data.  The  tests  were  con¬ 
ducted  for  periods  of  time  ranging  from  five-minute  segments  at  the  beginning  of 
a  run  to  20-minute  segments  after  one  hour  of  exposure. 

c.  Test  Sound  Pressure  Levels.  Test  sound  pressure  levels  were  chosen  such  that  the 
computed  minimum  time- to -failure  would  be  approximately  one  hour  for  the 
shortest  time  and  ten  hours  the  maximum  time.  The  upper  time  limit  was  varied 
when  necessary  to  keep  non-linearities  at  a  minimum. 

The  spectra  were  shaped  for  fatigue  tests  to  concentrate  the  acoustical  energy  in 
a  desired  range  of  frequencies*  It  was  determined  from  Reference  10  that  the 
bandwidth  for  testing  should  be  a  minimum  of  three  times  the  bandwidth  of 
response  in  order  for  the  structure  to  be  excited  with  broad -band  noise  relative  to 
the  response  bandwidth. 

d.  Test  Results.  Results  of  the  ambient  temperature  tests  will  be  described  beginning 
with  mode  studies  and  progressing  systematically  through  the  fatigue  tests. 

•  Mode  Shapes.  Mode  shapes  (Chladni  patterns)  for  each  panel  configuration 
were  determined  through  the  (5,5)  mode  for  both  the  0.051  and  0.032-inch 
thick  test  specimens.  Figures  8  through  10  are  photographs  of  the  mode  shapes 
for  an  0,051  M  panel  with  48-inch  radius,  an  0.03211  panel  with  48-Inch 
radius,  and  an  Q,032M  flat  panel.  Figure  8  should  be  compared  to  Figure  10 
to  get  the  effect  of  panel  thickness  on  frequency.  Figure  9  should  be  com¬ 
pared  to  Figure  10  to  get  the  effect  of  radius  of  curvature  on  frequency. 

The  mode  shape  studies  were  also  used  to  obtain  a  rigidly  controlled  check  on 
the  frequency  theory  derived  in  Section  II. A* 

•  Frequency  Sweeps.  Frequency  sweeps  were  made  at  various  sound  pressure 
levels  to  determine  response  characteristics,  damping  ratio,  and  range  of 
linear  response.  Figures  11  through  14  are  typical  plots  of  strain  response  for 
sinusoidal  frequency  sweep  excitation  at  130  decibels.  Frequency  was  swept 
from  50  cps  to  2000  cps  but  only  the  50  cps  to  600  cps  range  is  shown  because 
the  response  of  each  test  specimen  was  negligible  at  frequencies  above  600  cps. 
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FIGURE  6.  TYPICAL  STRAIN  GAGE  LOCATIONS  FOR  DETERMINING 
PANEL  RESPONSE 
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FIGURE  7.  SIMPLE  PANEL  TEST  FIXTURE  AND  TEST  SPECIMENS  INSTALLED  FOR  TEST 
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FIGURE  8.  MODE  SHAPES  FOR  0.051  PANEL  WITH  48-INCH  RADIUS 
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FIGURE  9.  MODE  SHAPES  FOR  0.032  FLAT  PANEL 
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FIGURE  10.  MODE  SHAPES  FOR  0.032  PANEL  WITH  48-INCH  RADIUS 
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FIGURE  11.  SINE  SWEEP  EXCITATION  -  FLAT  PANEL 
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FIGURE  12.  SINE  SWEEP  EXCITATION  -  96 
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FIGURE  13  SINE  SWEEP  EXCITATION-72' 
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FIGURE  14.  SINE  SWEEP  EXCITATION -48 


The  damping  ratio,  6,  was  determined  for  each  test  specimen  configuration 
using  the  bandwidth  method  and  the  log-decrement  method.  Both  methods 
produced  damping  ratios  which  compared  reasonably  well.  The  values  of  the 
damping  ratios  ranged  from  1%  to  4%.  These  lie  within  the  range  of  damping 
ratios  for  simple  panels  published  in  Reference  1 1  . 

•  Broad-Band  Excitation.  Each  configuration  of  curved  panel  was  exposed  to 
broad-band  noise  having  a  flat  spectrum  from  50  to  2,000  cps.  Dynamic 
strain  measurements  were  made  of  the  most  active  gage  (gage  ^1)  for  each 
configuration.  The  strain  data  were  recorded  on  magnetic  tape  and  later 
analyzed  by  a  narrow-band  (1  cps)  analyzer  to  compare  the  sinusoidal 
frequency  sweep  responses  with  the  flat  spectrum,  broad-band  excitation 
responses.  Figures  15  through  18  are  nominal  one-eps  bandwidth  analyses  of 
strain  responses  resulting  from  flat  spectrum,  broad-band  noise  excitation  at 
an  overall  sound  pressure  level  of  151  decibels. 

•  Brood-Band  Acoustical  Noise  Spectra,  In  a  previous  section  it  was  pointed 
out  that  a  shaped  spectrum  can  be  used  for  broad-band  testing  provided  the 
bandwidth  of  the  excitation  is  at  least  three  times  the  bandwidth  of  the  response. 
The  acoustical  noise  spectra  for  the  fatigue  tests  were  shaped  accordingly  and 
are  shown  in  Figure  19.  Two  spectra  were  required  to  satisfactorily  encompass 
the  significant  response  peaks  (l ,  1  mode)  of  the  four  panel  configurations. 

The  spectrum  of  Figure  19  drawn  with  a  dashed  line  was  used  to  test  the  flat 
panels  and  those  with  the  48-inch  radius.  The  spectrum  drawn  with  a  solid 
line  was  used  to  test  the  panels  with  the  72-inch  radius  and  the  96-inch 
radius.  These  spectra  were  raised  or  lowered  to  obtain  a  desired  test  level . 

There  was  no  noticeable  change  in  spectrum  shape  when  the  test  levels  were 
raised  or  lowered  as  required  for  the  fatigue  tests. 

•  Fatigue  Tests.  Twenty-four  test  specimens  of  the  designs  shown  in  Figure  4 
were  tested  to  failure  at  the  test  spectrum  levels  shown  In  the  following  table. 


TABLE  4 


SIMPLE  SPECIMEN  TEST  SUMMARY 


Configuration 

Replicates 

Test  Spectrum  Level, 
db 

(re.  0.0002  Microbar) 

Flat 

2 

134 

Flat 

2 

137 

Flat 

2 

140 

96"  R 

2 

136 

96"  R 

2 

139 

96"  R 

2 

142 

72"  R 

2 

135.5 

72"R 

2 

138.5 

72"R 

2 

141.5 

48"  R 

2 

136 

48"  R 

2 

139 

48"  R 

2 

142 

34 


CO 

o_ 

O 


> 

U 

z 

LU 

ID 

o 

LU 

Oc- 


1VN9IS  HV3  JO  !N3Dil3d-NIVyiS  SWH 


35 


FIGURE  15.  BROAD-BAND  NOISE  EXCITATION,  151  DB  -  FLAT  PANEL 
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FIGURE  16.  BROAD-BAND  NOISE  EXCITATION,  151  DB  -  96' 
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FIGURE  17.  BROAD-BAND  NOISE  EXCITATION,  151  DB  -  72' 
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FIGURE  18.  BROAD-BAND  NOISE  EXCITATION,  151  DB  -  48 
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FIGURE  19.  BROAD-BAND  NOISE  SPECTRA  FOR  FATIGUE  TESTS 


The  results  of  the  fatigue  tests  are  presented  in  the  form  of  fatigue  curves, 
Figure  20.  The  fatigue  curves  were  plotted  using  spectrum  level,  db,  as  the 
ordinate  and  mean  cycles-to-failure  as  the  abscissa,  (Cycles- to-fai lure  were 
determined  as  the  product  of  time-to-failure  and  critical  response  frequency.) 

Figure  21  is  presented  to  show  the  relationship  between  maximum  root-mean- 
square  (rms)  panel  strain  versus  eye ies-to-fai lure .  As  expected,  these  e  -  N 
curves  have  the  same  general  appearance  as  the  curves  of  Figure  20. 

Magnetic  tape  recordings  of  strain  which  produced  response  spectra  of  the 
type  shown  in  Figures  15  through  18  were  analyzed  to  determine  the  nature  of 
their  statistical  properties.  A  30-second  sample  of  the  maximum  strain 
response  from  each  of  the  four  test  specimen  configurations  was  analyzed  by  a 
probability  analyzer  to  determine  probability  distribution  of  strain  peaks.  The 
strain  data  analyzed  were  those  obtained  at  the  highest  test  spectrum  level. 
Figure  22  is  a  plot  of  the  probability  distribution  of  the  strain  peaks  for  the 
four  configurations. 


4.  Elevated  Temperature  Tests 

The  object  of  this  phase  of  the  simple  specimen  experimental  investigation  was  to  determine 
the  relative  fatigue  damage  resulting  from  alternate  application  of  heat  and  acoustical 
excitation  and  simultaneous  application  of  heat  and  acoustical  excitation. 

A  description  of  the  test  specimens,  test  setup,  test  procedure,  and  test  results  will  follow 
In  that  order. 

a.  Test  Specimen.  The  test  specimens  were  flat  with  overall  dimensions  of  11  .0"  x 
13.0"  x  0.032"  .  Four  of  the  specimens  were  8AI-IM0-IV  titanium  and  four  were 
17-7  PH  stainless  steel.  Holes  for  3/16-inch  diameter  bolts  were  drilled  around 
the  perimeter  of  each  specimen  as  shown  on  the  drawing  of  the  test  specimen. 
Figure  4. 

b.  Test  Setup.  Each  test  specimen  was  attached  to  a  picture  frame  mounting  fixture 
of  the  type  shown  in  Figure  5.  A  torque  of  45-55  inch-pounds  was  applied  to 
each  bolt  to  keep  the  edge  fixity  uniform  for  all  specimens. 

Two  of  the  mounting  fixtures  were  bolted  to  each  side  of  the  progressive  wave 
test  section  of  the  High  Intensity  Sound  System  as  shown  in  Figure  23.  This  photograph 
shows  that  a  removable  side  of  the  test  section  has  been  opened  for  access.  Also, 
three  monitor  microphones  are  located  along  the  vertical  centerline  of  the  test 
specimens  and  one  located  forward  at  the  mouth  of  the  horn. 

Quartz  lamp  tubes  were  used  to  supply  the  heat  for  the  test  specimens.  Figure  24 
is  a  photograph  of  these  heaters  moved  away  from  two  test  specimens  to  facilitate 
inspection  for  fatigue  crocks  and  permanent  thermal  buckling. 

Each  test  specimen  was  instrumented  with  thermocouples  to  determine  the 
temperature  distribution  during  elevated  temperature  tests.  Test  specimen 
thermocouple  locations  are  shown  in  Figure  25.  The  thermocouples  were  held  to 
the  surface  of  each  test  specimen  by  aluminized  tape. 

High  temperature  strain  gages  were  mounted  on  the  test  specimens  In  the  pattern 
shown  in  Figure  26. 
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FIGURE  20.  SPECTRUM  LEVEL  VERSUS  MEAN  CYC LES-TO-FAI LURE 
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FIGURE  21  .  €  -N  CURVES  FOR  SIMPLE  CURVED  PANELS 


PROBABILITY  OF  PEAKS  EXCEEDING 


FIGURE  22.  PROBABILITY  DISTRIBUTION  OF  STRAIN  PEAKS 
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FIGURE  23.  VIEW  OF  TEST  SPECIMENS  MOUNTED  IN  TEST  SECTION 
(A  REMOVABLE  SIDE  HAS  BEEN  OPENED  FOR  ACCESS) 
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FIGURE  24.  QUARTZ  LAMP  HEATERS  FOR  ELEVATED  TEMPERATURE  TESTS 
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FIGURE  25.  THERMOCOUPLE  LOCATIONS 
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17-7  PH  STAINLESS  STEEL 


1.75 

* 

1.75 

_ A _ 

°T 

2  T 

1.75-*| 

1  0 

3ca 

SI,  S2,  S3 

S4 

8AI-IM0-IV  TITANIUM 


1  .0 

_ 4 _ 

4Q  —f 

1 

3f 

2ea 

T1 


FIGURE  26.  STRAIN  GAGE  LOCATIONS  -  ELEVATED  TEMPERATURE  TESTS 
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c*  Test  Procedure,  Frequency  sweeps  were  made  at  a  number  of  discrete  frequency 
sound  pressure  levels  to  determine  the  response  characteristics  of  the  two  con¬ 
figurations.  These  frequency  sweeps  were  made  at  ambient  temperature  and  test 
temperature  of  450°  F . 

The  frequency  sweeps  were  used  to  shape  broad-band  sound  pressure  test  spectra 
for  the  fatigue  tests.  Results  of  the  ambient  temperature  tests  reported  in  Section 
II. B. 3  showed  similar  test  specimen  response  characteristics  for  broad-band  exci¬ 
tation  and  discrete  frequency  excitation.  Consequently,  the  low-intensity, 
broad-band  tests  were  omitted  for  this  phase  of  the  simple  test  specimen  investi¬ 
gation  . 

The  procedure  for  the  fatigue  test  on  one  set  of  four  specimens  was  as  follows: 

Two  specimens,  say  1  and  2,  were  exposed  to  heat  generated  by  quartz  lamps. 
Shortly  after  reaching  thermal  equilibrium,  the  sound  generators  were  started  and 
the  sound  pressure  level  brought  up  to  the  test  level  on  all  four  specimens.  The 
two  specimens,  3  and  4,  mounted  on  the  opposite  side  of  the  test  section  were 
excited  acoustically  at  ambient  temperature.  After  a  predetermined  time  of 
exposure  the  acoustical  excitation  and  heat  were  removed.  The  heaters  were  then 
moved  to  specimens  3  and  4.  They  were  heated  to  the  same  temperature  and  for 
the  same  length  of  time  as  specimens  1  and  2,  but  without  acoustical  excitation. 
This  completed  one  cycle  of  the  fatigue  test.  The  specimens  were  dye-penetrant 
checked  for  fatigue  cracks  after  each  acoustical  excitation.  When  a  crack  was 
discovered  its  location  and  length  were  recorded  versus  test  time  until  the  panels 
were  removed  fnom  the  test  chamber. 

d*  Test  Results*  Results  of  the  elevated  temperature  tests  will  be  described.  The 
description  will  begin  with  the  frequency  sweeps  and  end  with  the  fatigue  tests  of 
the  titanium  and  stainless  steel  specimens. 

•  Frequency  Sweeps,  Figures  27  through  29  are  plots  of  strain  response  for 
sinusoidal  frequency  sweep  excitation  at  130  decibels.  Frequency  sweeps  for 
the  titanium  and  stainless  steel  specimen  were  made  of  450°  F  and  ambient 
temperature. 

Damping  was  determined  for  each  panel  configuration  and  test  condition. 
Damping  ratios  determined  by  the  bandwidth  method  range  from  1 .3%  to 
2.1%.  The  damping  ratios  at  elevated  temperatures  were  less  than  ambient 
temperature  ratios. 

•  Broad-Band  Acoustical  Noise  Spectra.  It  was  necessary  to  shape  the  test 
spectrum  so  that  fatigue  failures  could  be  produced  in  a  reasonable  length  of 
time.  Figure  30(b)  is  a  typical  test  spectrum  for  the  elevated  and  room  temper¬ 
ature  tests  of  the  titanium  specimens.  Figure  30(a)  is  a  typical  test  spectrum 
for  the  elevated  and  room  temperature  tests  of  the  stainless  steel  specimens. 

The  spectrum  for  the  stainless  steel  specimens  had  two  lobes.  One  lobe  was 
required  to  excite  specimen  response  at  ambient  temperature  whereas  the  other 
lobe  was  needed  to  excite  specimen  response  at  elevated  temperature. 

•  Temperature  ,  A  table  of  typical  temperature  distribution  over  the  specimen 
surfaces  during  frequency  sweeps  and  fatigue  tests  Is  shown  below. 
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FIGURE  27.  SINE  SWEEP  EXCITATION,  AMBIENT  TEMPERATURE-TITANIUM 
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FIGURE  28.  SINE  SWEEP  EXCITATION,  450°  F  -  TITANIUM 
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FIGURE  29.  SINE  SWEEP  EXCITATION,  AMBIENT  AND  450°  F  -  STAINLESS  STEEL 
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FIGURE  30.  BROAD-BAND  NOISE  SPECTRA  FOR  ELEVATED  TEMPERATURE  TESTS 


TABLE  5 

TYPICAL  FAHRENHEIT  TEMPERATURE  DISTRIBUTION 


Thermocouple 

No . 

(See  Figure  25) 

Ti-8AI-lMo-lV 

17-7PH  Stainless  Steel 

Frequency 

Sweeps 

Fatigue 

Tests 

Frequency 

Sweeps 

Fatigue 

Tests 

1 

430 

420 

410 

460 

2 

420 

460 

420 

470 

3 

490 

460 

490 

420 

4 

420 

390 

460 

460 

5 

450 

470 

470 

520 

6 

430 

430 

- 

- 

7 

450 

460 

- 

480 

8 

- 

450 

- 

440 

9 

470 

440 

- 

390 

10 

420 

420 

420 

490 

11 

480 

470 

420 

470 

12 

440 

490 

- 

— 

•  Fatigue  Tests.  Four  titanium  and  four  stainless  steel  test  specimens  were 
tested  to  failure  at  a  spectrum  level  of  141  db  using  the  test  procedure 
described  previously.  The  results  of  these  tests  are  contained  in  Tables  6  and 
7.  Sketches  of  crack  location  and  tabulated  crack  propagation  are  also 
included . 


C.  Comparison  of  Theoretical  and  Experimental  Results 


Natural  frequencies  for  one  of  the  test  panels  described  In  Section  II. B  were  determined 
experimentally  to  assess  the  accuracy  of  the  theoretical  frequency  analysis.  The  edge 
conditions  of  the  test  setup  described  in  Section  II. B  could  be  classified  as  being  some¬ 
where  between  clamped  and  simply-supported.  In  an  effort  to  make  the  edges  consistent 
with  theory  the  frame  shown  in  Figure  5  was  reinforced  with  a  1/2"  thick  aluminum  picture 
frame  type  member  bolted  to  the  outer  surface  of  the  panel .  Extra  fasteners  were  added  to 
increase  the  clamping  effects. 

The  mode  shapes  and  frequencies  were  determined  as  follows:  First,  two  loudspeakers 
which  could  be  freely  positioned,  were  used  to  acoustically  excite  the  different  modes. 

One  of  the  speakers  was  equipped  with  a  phase  reversal  switch  in  the  line  because  even 
modes  were  most  easily  excited  with  the  speakers  out-of-phase. 

Table  8  shows  calculated  values  of  natural  frequency  for  both  clamped  and  simply -supported 
edge  conditions  along  with  the  experimental  values  for  a  9"  x  11"  aluminum  panel  having  a 
96"  radius  of  curvature  and  0.048"  thickness.  The  comparison  is  quite  good  except  for  the 
(1,1)  through  (1,4)  modes.  This  difference  probably  lies  in  the  matching  of  boundary  con¬ 
ditions  between  theory  and  experiment. 
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TABLE  6 

FAILURE  SUMMARY:  TITANIUM  SPECIMENS 


HEAT  AND  SOUND  SIMULTANEOUSLY 


TEST 

TIME 

MIN, 

CRACK 

NUMBER 

CRACK 

LENGTH 

INCHES 

120 

1 

1.3 

— 

— 

— 

— 

— 

— 

TEST 

TIME 

MIN. 

CRACK 

NUMBER 

C 

LI 

| 

;rac 

ENG 

K 

TH 

Eft 

20  1 

1 

t 

S 

1.8 

0.6 

30 

1 

2 

3 

25 

OS 

0.6 

45 

1 

2 

3 

3,0 

3.0 

3 £> 

*  ■  *  '  A. . 

2 

•  * 

•  * 

l  • 

’  1  v 

T  3  3> 

• 

•  ■ 

V  * 

: . L---: 

t  4  ; 

1  1 

•  0 

•  1 

HEAT  AND  SOUND  ALTERNATELY 


TEST 

TIME 

MIN. 

CRACK 

NUMBER 

CRACK 

LENGTH 

INCHES 

420 

1 

2 

3 

0.1 

02 

,  * 

540 

1 

2 

3 

- 

02 

— 

645 

1 

t 

3 

3  JO 

0.3 

1.0 

TEST 

CRACK 

CRACK 

TIME 

NUMBER 

LENGTH 

MIN. 

INCHES 

645 

1 

0.1 

— 

— 

— 

— 

— 

- - 
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TABLE  7 

FAILURE  SUMMARY:  STAINLESS  STEEL  SPECIMENS 


!  si  : 

(  * 

*  1 

*  , 

;  s  2  ; 

§  ( 

*  i 

• 

i 

HEAT  AND  SOUND  SIMULTANEOUSLY 


TEST 

TIME 

MIN, 

CRACK 

NUMBER 

CRACK 

LENGTH 

INCHES 

30 

t 

0.03 

60 

1 

0,75 

90 

1 

1.75 

TEST 

TIME 

MIN. 

CRACK 

NUMBER 

CRACK 

LENGTH 

INCHES 

90 

1 

1.50 

— 

— 

— 

— 

— 

— 

HEAT  AND  SOUND  ALTERNATELY 


TEST 

TIME 

MIN. 

CRACK 

NUMBER 

CRACK 

LENGTH 

INCHES 

200 

1 

3.0 

200 

2 

0.75 

200 

3 

0.10 

TEST 

CRACK 

CRACK 

TIME 

NUMBER 

LENGTH 

MIN. 

INCHES 

- NO  1 

FAILURES 

TEST  STOPPED  AFTER 
SSO  MIN.  TEST  TIME 
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TABLE  8 

FREQUENCY  STUDIES  ON  CURVED  PANEL  WITH  96"  RADIUS 


n 

m 

1 

2 

3 

4 

5 

346 

471 

778 

1225 

1795 

Clamped -Theory 

i 

231 

385 

696 

1128 

1746 

Experimental 

Simply  Supported  Theory 

165 

273 

557 

960 

1478 

440 

590 

900 

1347 

1917 

Clamped-Theory 

2 

436 

561 

861 

1283 

1815 

Experimental 

323 

407 

677 

1077 

1594 

Simply  Supported  Theory 

628 

794 

1104 

1548 

2117 

Clamped-Theory 

3 

636 

793 

1077 

1492 

2025 

Experimental 

496 

611 

876 

1272 

1788 

Simply  Supported  Theory 

909 

1085 

1393 

1832 

2397 

Clamped-Theory 

4 

899 

1093 

1371 

1775 

2297 

Expe  rimental 

741 

882 

1149 

1544 

2059 

Simply  Supported  Theory 

1278 

1458 

1765 

2199 

2759 

Clamped-Theory 

5 

1245 

1441 

1728 

2125 

— 

Experimental 

1069 

1224 

1498 

1893 

2407 

Simply  Supported  Theory 

D  .  Discussion  of  Results 


1  .  Non-Dimensional  Frequency  Charts 

The  frequency  analysis  for  isotropic  curved  panels  with  no  coupled  mode  terms  (equation 
(30))  has  been  programmed  in  Fortran  language  for  solution  on  the  IBM  7094,  The  equa¬ 
tions  are  non-dimensionalized  in  terms  of  three  independent  variables.  A,  cp,  and  t,  and 
the  dependent  variable  which  is  non-dimensional  frequency. 

Based  upon  information  related  to  practical  aircraft  structure,  a  range  of  the  aforementioned 
variables  was  selected  as  follows: 


0  £  b/a  =  cp  ^  3.14 

20  ~  t  s  1000 

0.5  s  b/l  =  A  s  2.0 

Values,  including  the  limits,  of  each  of  the  above  parameters  were  selected  and  frequency 
was  calculated  for  all  combinations  for  both  simply  supported  and  clamped  edges. 

For  particular  values  of  aspect  ratio,  A,  non-dimensional  frequency  is  plotted  for  six 
modes  and  six  values  ofcp  for  variation  of  length-to-thickness  ratio .  Figures  31  to  35  are 
for  simply  supported  edges  and  Figures  36  to  40  give  clamped  edge  frequencies. 
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FIGURE  31 .  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
SIMPLY-SUPPORTED  EDGES,  A  =  0.50 
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FIGURE  32.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
SIMPLY-SUPPORTED  EDGES,  A  =  0.67 
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FIGURE  33.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
SIMPLY-SUPPORTED  EDGES,  A  =  1  .00 
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FIGURE  34.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
SIMPLY-SUPPORTED  EDGES,  A  =  1.50 
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FIGURE  35.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
SIMPLY-SUPPORTED  EDGES,  A  =  2.00 
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FIGURE  36.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
CLAMPED  EDGES,  A  =  0.50 


62 


1000 


400 


100 


40 


10 


Code 

Mode  Number 

Straight  Curved 
Edge,  m  Edge,  n 

Code 

Subtended 
Angle,  4 

c 

1 

1 

1 

0 

d 

1 

2 

2 

.05 

e 

1 

3 

3 

.1 

f 

2 

1 

4 

.2 

a 

2 

2 

5 

1.0 

B 

h 

3 

1 

6 

3.14 

J  I  L 


J — L 


10 


40  100 

length  to  thickness  ratio,  — 

h 


400 


1000 


FIGURE  37.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
CLAMPED  EDGES,  A  =  0.67 
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FIGURE  38.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
CLAMPED  EDGES,  A  =  1  .00 
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FIGURE  39.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
CLAMPED  EDGES,  A  =  1  .50 
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FIGURE  40.  NON-DIMENSIONAL  FREQUENCY  SOLUTIONS 
CLAMPED  EDGES,  A  =  2.00 
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Once  non-dimensional  frequency  is  found  then  actual  frequency  can  be  determined  from 
the  nomogram  (Figure  41). 

For  example,  calculate  the  natural  frequencies  of  a  clamped,  curved  panel  having  the 
following  dimensions: 

radius,  a  —  1 0QM 

arc  length,  b  =  10" 
length,  t  =  20" 

thickness,  h  =  0.05" 


The  non-dimensional  ratios  are: 


A  =  0.5 

cp  =  0.1 

t  =400 

Table  9  shows  values  of  X.  for  the  different  combinations  of  mode  number.  These  values 
were  taken  from  Figure  36  for  A  =  0.5.  Frequency  converted  through  the  use  of  the  nomo¬ 
gram  is  also  displayed  in  Table  9. 


TABLE  9 

NATURAL  FREQUENCIES  FOR  EXAMPLE  PROBLEM 


m 

n 

\ 

f 

i 

i 

51 

300 

i 

2 

65 

382 

i 

3 

101 

594 

2 

1 

54 

318 

2 

2 

71 

418 

3 

1 

61 

359 

2.  Experimental 

a .  Ambient  Temperature  Tests. 

•  Mode  Shapes.  Some  of  the  modes  could  not  be  excited .  From  Figure  8  it  can 
be  seen  that  the  (4,2),  (4,4),  and  (5,4)  modes  were  unresponsive.  Also,  the 
(5,5)  mode  was  unresponsive  for  the  0. 032-Inch  fiat  panel  as  shown  in  Figure 
9.  Figure  10  shows  tnat  seven  modes,  (2,1),  (2,3),  (2,4),  (4,4),  (4,5),  ^(5,3), 
and  (5,4)  did  not  readily  respond.  The  lack  of  response  is  attributed  to  high 
damping  and  low  excitation. 
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FIGURE  41  .  NOMOGRAM  FOR  CONVERTING  NON-DIMENSIONAL  FREQUENCY  TO  ACTUAL  FREQUENCY 


Cross-over  effects,  (1,2)  mode  having  lower  frequency  than  the  (1,1)  mode, 
were  observed  for  each  of  the  curved  panels  with  0. 032-inch  thickness, 

•  Frequency  Sweeps.  Strain  activity  above  the  first  mode  of  response  was 
negligible  for  alt  specimen  configurations  except  the  panel  with  the  48-inch 
radius.  It  showed  some  activity  at  665  cps  which  is  the  (1 ,3)  mode.  The 
strain  associated  with  this  higher  mode  was  approximately  30  percent  of  the 
strain  at  the  fundamental  frequency. 

Unlike  the  mode  shape  studies,  cross-over  effects  were  not  observed  during 
the  frequency  sweeps.  Figure  12  through  14.  Out-of-phase  speakers  were 
required  to  excite  the  (1,2)  mode  during  the  mode  shape  studies.  Therefore, 
it  is  possible  that  the  acoustical  excitation  and  test  specimen  did  not  have  the 
required  phase  relationship  to  excite  the  (1,2)  mode  when  the  sinusoidal 
frequency  sweeps  were  made. 

Significant  frequencies  determined  from  the  mode  shape  studies  did  not 
compare  with  those  determined  from  the  frequency  sweeps.  This  lack  of 
correlation  is  attributed  to  dissimilar  edge  fixities  and  ambient  temperatures. 

Sub-harmonics  of  the  type  described  by  Rucker  in  Reference  12  were  not 
observed  at  any  time. 

•  Broad-Band  Excitation.  Strain  spectra  resulting  from  broad-band  acoustical 
excitation.  Figures  15  through  18,  had  the  same  general  shape  as  the 
sinusoidal  frequency  sweeps  in  the  vicinity  of  the  (1,1)  mode,  but  were  more 
active  at  the  higher  frequencies.  This  activity  in  the  higher  modes  was  con¬ 
sidered  insignificant  relative  to  fatigue  damage  because,  with  the  exception 
being  the  72"  radius,  the  higher  mode  strain  magnitudes  were  less  than  30 
percent  of  the  strain  response  of  the  (1,1)  mode . 

Cross-over  effects  were  observed  when  the  curved  panels  were  exposed  to 
broad-band  excitation.  Evidence  of  this  effect  is  present  in  Figures  16 
through  18  where,  in  each  case,  a  small  crinkle  in  the  response  curve  pre¬ 
cedes  the  maximum  response  peak.  This  crinkle  Is  indicative  of  the  (1,2) 
mode  response. 

•  Fatigue  Tests.  It  can  be  seen  from  Figure  21  that  curvature  hod  a  pronounced 
effect  on  fatigue  resistance  with  the  degree  depending  on  the  radius  of  curva¬ 
ture.  The  test  specimens  with  a  96-inch  radius  of  curvature  showed  only  a 
slight  advantage  over  the  flat  test  specimens.  The  test  specimens  with  the  72- 
inch  radius  and  the  48-inch  radius  demonstrated  considerably  more  fatigue 
resistance  than  the  flat  test  specimens.  At  a  spectrum  level  of  139  db,  the 
average  eye !es-to-fai lure  advantage  of  the  test  specimens  having  radii  of 
curvature  of  96  inches,  72  inches,  and  48  inches  over  the  flat  test  specimen 
was  determined  from  Figure  20  to  be  1 .2,  1 .7,  and  5  times  respectively. 

The  slopes  of  the  fatigue  curves.  Figure  20,  for  the  test  specimens  having 
radii  of  infinity,  96,  and  72  inches  were  the  same  for  practical  purposes.  The 
flatter  slope  for  the  test  specimen  with  the  48-inch  radius  can  be  attributed  to 
data  scatter. 

The  e  -  N  curves  of  Figure  21  show  that  tests  specimen  fatigue  resistance 
increases  inversely  as  the  radius  of  curvature. 

Cross-plots  of  rms  versus  1/R,  R  being  the  radius  of  curvature  in  feet,  and  1/R 
versus  mean  eye les-to-fai lure  are  shown  in  Figures  42  and  43.  The  curves  of 
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FIGURE  42.  1 /R  VERSUS  RMS  STRAIN 
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CYCLES  -  N 

FIGURE  43.  1/R  VERSUS  MEAN  CYC LES-TO-FAI LURE 


Figure  43,  fatigue  life  of  the  test  specimens  as  a  function  of  curvature, 
clearly  show  that  the  most  significant  increase  in  fatigue  life  is  the  increase 
of  radius  from  72  inches  to  48  inches.  This  does  not  agree  with  the  conclu- 
sions  drawn  from  discrete  frequency  tests  of  curved  panels  reported  in  Reference 
13.  A  possible  explanation  of  the  difference  in  results  of  Reference  13  and 
this  investigation  is  the  nature  of  the  failures  experienced  in  each  case.  In 
Reference  13,  it  was  stated  that  the  flat  panels  failed  near  the  bolt  heads 
whereas  the  curved  panels  generally  failed  along  the  long  side  near  the  edge 
of  the  supporting  frame.  In  this  Investigation,  the  nature  of  failure  for  both 
the  flat  panels  and  the  curved  panels  was  the  same.  Fatigue  cracks  formed  in 
every  case  along  the  fastener  row  as  shown  in  Figure  44.  The  supporting 
frame  had  rounded  edges  to  prevent  cracks  from  forming  there. 

The  experimentally  determined  probability  distribution  of  strain  peaks^for  each 
pane!  configuration  as  shown  in  Figure  22  lies  very  close  to  the  Rayleigh  dis¬ 
tribution  line  for  peak-to-rms  ratios  up  to  three.  The  probability  distribution 
of  strain  peaks  will  deviate  appreciably  from  the  Rayleigh  distribution  if  the 
response  is  nonlinear  and  multiple-modal.  Therefore,  data  plotted  on  Figure 
22  show  that  the  response  of  each  test  specimen  configuration  was  essentially 
unimodal  and  linear. 

•  Strain  Measurements.  Response  of  the  test  specimens  to  the  acoustical 
excitation  was  determined  through  the  use  of  electrical  strain  gages  attached 
to  the  specimen  at  points  of  high  strain.  The  strain  magnitudes  thus  measured 
can  only  be  used  qualitatively  because  the  placement  and  orientation  of  a 
strain  gage  permits  sensing  the  strain  only  at  one  small  area  and  only  In  the 
direction  in  which  the  gage  is  sensitive.  Only.if  the  gage  is  placed  where 
failure  ultimately  occurs  and  if  the  gage  direction  Is  the  principal  strain^ 
direction  for  the  dominant  mode  at  this  point  is  the  measured  strain  magnitude 
quantitatively  meaningful.  Therefore,  strain  magnitudes  shown  in  this  report 
are  nominal  values  only. 

b.  Elevated  Temperature  Tests. 

•  Frequency  Sweeps.  Sinusoidal  frequency  sweeps  for  the  titanium  test  specimens 
at  ambient  and  elevated  temperatures.  Figure  27  and  28,  show  negligible 
response  above  the  (1,1)  mode.  The  frequency  of  the  (1,1)  mode  at  ambient 
temperature  was  140  cps  whereas  the  frequency  of  this  mode  shifted  downward 
to  95  cps  at  450°  F.  Also,  the  strain  magnitude  at  450°  F  was  60  percent 
greater  than  the  strain  response  at  ambient  temperatures.  It  is  believed  that 
the  heated  titanium  panel  strain  response  was  strongly  influenced  by  the  first 
thermal  buckling  mode*  This  mode  Is  characterized  by  a  general  bowing  of 
the  panel  in  one  direction.  When  heat  was  applied,  the  titanium  panels  were 
observed  to  have  a  buckling  pattern  similar  to  that  of  the  first  thermal  buckling 
mode  * 

Frequency  sweeps  for  the  stainless  steel  panels  were  similar  to  those  for 
titanium.  At  ambient  temperature  there  was  negligible  activity  above  the 
(1,1)  mode.  Figure  29,  which  occurred  at  130  cps.  At  450  F  the  maximum 
response  appeared  at  265  cps  as  shown  in  Figure  29.  The  strain  at  this 
frequency  was  two  times  greater  than  the  maximum  strain  response  ot  ambient 
temperature.  There  was  some  response  at  385  cps  but  it  was  considered  insig¬ 
nificant  relative  to  its  contribution  to  fatigue  damage. 
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FIGURE  44.  TYPICAL  FATIGUE  CRACK  FOR  CURVED  TEST  SPECIMENS 


The  elevated  temperature  response  of  the  stainless  steel  test  specimen  was  not 
the  same  mode  as  that  of  the  titanium.  The  maximum  response  frequency 
shifted  upward  from  130  cps  to  265  cps.  This  upward  shift  in  frequency  was 
probably  caused  by  second-mode  thermal  buckling  rather  than  first-mode 
buckling . 

•  Temperatures.  Considerable  difficulty  was  encountered  in  obtaining  uniform 
temperature  over  the  surface  of  the  test  specimens.  This  was  due  primarily  to 
three  things:  The  area  of  the  test  specimens  was  small,  the  mounting  frame  was 
a  heat  sink  in  spite  of  the  thin  phenolic  spacer  placed  between  the  frame  and 
the  test  specimen,  and  the  air  flow  in  the  test  section,  necessary  in  the  sound 
generation,  caused  uneven  convection  heat  losses  over  the  test  specimen 
surface.  A  one-inch  width  of  heat-resistant  black  paint  was  applied  to  the 
edges  of  the  test  specimen  to  increase  the  temperature  around  the  periphery 
where  the  heat  loss  was  the  greatest. 

•  Fatigue  Tests.  A  comparison  of  the  time  to  first  failure  for  the  specimens 
exposed  to  heat  and  sound  simultaneously  to  that  of  the  specimens  exposed  to 
heat  and  sound  alternately  revealed  that  the  former  was  considerably  more 
severe  than  the  latter.  This  is  true  for  both  titanium  and  stainless  steel. 

Based  upon  the  data  contained  in  Tables  6  and  7,  the  fatigue  life  of  the 
titanium  specimens  exposed  to  heat  and  sound  simultaneously  was  about  one- 
eighth  the  fatigue  life  of  the  specimens  exposed  to  heat  and  sound  alternately. 
The  ratio  for  the  stainless  steel  specimens  was  about  one-fourth. 

Crack  propagation  was  also  greater  for  the  combined  environmental  condition. 

Figures  45a  and  45b  are  typical  strain  spectra  for  the  titanium  and  stainless 
steel  specimens  during  the  combined  heat  and  acoustical  excitation  tests. 

Figure  45a  shows  the  titanium  panel  responded  in  the  first  thermal  buckling 
mode  with  very  little  activity  at  higher  frequencies.  The  stainless  steel  panel 
experienced  its  maximum  response  at  the  second  thermal  buckling  mode, 

Figure  45b,  with  insignificant  response  at  other  frequencies. 
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FIGURE  45.  BROAD-BAND  NOISE  EXCITATION,  162  DB  -  450°  F 


Ill 

COMPLEX  SPECIMEN  INVESTIGATION 
A.  Theoretical  and  Analytical 


1  ,  Introduction 

Honeycomb  sandwich  structure  has  been  used  in  the  construction  of  aerospace  vehicles  for 
the  past  decade.  This  type  of  construction  is  generally  used  because  environmental  prob¬ 
lems  will  not  permit  usage  of  conventional  aircraft  structure  due  to  weight  penalties. 
Honeycomb  sandwich  is  the  best  structure,  weight-wise,  for  high  intensity  noise  environ¬ 
ments  because  of  its  resistance  to  acoustical  fatigue. 

Most  honeycomb  sandwich  structures  are  in  the  elemental  form  of  panels.  These  panels 
may  be  either  flat  or  curved  .  In  general,  the  design  of  these  panels  has  evolved  from  that 
of  a  uniform  thickness  to  the  more  efficient,  lighter  weight  structure  typical  of  that  shown 
in  Figure  49.  This  design  consists  of  a  core  bonded  to  two  thin  facing  sheets,  which  may 
be  of  different  thickness.  The  core  gradually  tapers  to  zero  thickness  near  the  edge  and 
the  inner  skin  stresses  are  transmitted  to  the  edge  through  a  bonded  doubler  and  "wet- 
layup"  fiberglass  edge  member. 

In  previous  vibration  analyses  of  honeycomb  sandwich  panels,  uniform  thickness  of  the  core 
and  face  sheets  has  been  assumed  (References  14  and  15).  In  these  analyses,  the  boundaries 
are  simply  supported . 

The  analysis  to  be  presented  in  this  section  includes  the  effects  of  the  tapered  edge  on 
natural  modes  of  vibration.  The  theoretical  model  chosen  to  represent  the  edge  is  shown 
in  Figure  46.  This  model  Is  compared  with  the  actual  edge  geometry.  Although  the  taper 
is  reduced  to  a  step-discontinuity  in  thickness,  it  is  felt  that  this  effect  on  the  theory  will 
be  negligible. 

The  frequency  analysis  which  follows  will  be  based  upon  a  Rayleigh-Ritz  energy  analysis. 
Mode  shapes  are  assumed  which  satisfy  the  geometric  boundary  conditions. 


2.  Honeycomb  Panel  Frequency  Analysis 
The  strain  energy  density  of  an  elastic  body  is 

U0  =  -g  LCTjJ  €e.}  (1) 

In  the  notation  introduced  by  Love  (Reference  2),  a,  =  a  ,  02  ”  o  ,  =  o^f  a ^  =  a  f 

a  =a  ,  a*  =a  -  The  notation  is  identical  forstrains. 

5  zx  6  xy 

The  Hooke's  law  relationship  between  stress  and  strain  is  then  expressed  as 

Lg.J  =  Le.J  C..  (2) 

u  f  u  r  tj 
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Theoretical  Model 
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FIGURE  46.  THEORETICAL  AND  ACTUAL  MODELS  OF  TAPERED  EDGE  GEOMETRY 
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For  an  orthotropic  linear  elastic  solid  the  elastic  coefficient  matrix,  [C.J,  is  symmetric 
and  has  the  following  form: 


'I 


C11 

C12 

C13 

0 

0 

0 

C21 

C22 

C23 

0 

0 

0 

C31 

C32 

C33 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

C55 

0 

0 

0 

0 

0 

0 

C12  =  C2J 
C13  =  C31 
C23  =  C32 


(3) 


Using  these  relationships,  the  strain  energy  density  becomes 

u0=iLeiJ  [pri] 


(4) 


In  expressing  strain  in  terms  of  displacement,  the  core  will  be  treated  as  a  thick  shell; 
however,  in  order  to  reduce  the  complexity  of  the  analysis  the  facing  sheets  will  be  con¬ 
sidered  as  thin  membranes. 

In  considering  the  facings  as  membranes,  it  is  implicitly  assumed  that  no  transverse  shear 
strain  occurs  (i  .e.,  cp  =  i|i  =  0).  To  make  this  assumption  since  there  are  minor  transverse 
shear  stresses,  the  elastic  coefficients  for  transverse  shear,  C44  and  C55,  are  assumed  to 
be  infinite.  Because  of  this  assumption,  C44  and  C55  must  be  omitted  from  the  elastic 
coefficient  matrix.  This  assumption  is  not  detrimental  to  the  analysis  since  facing  sheets 
are  usually  an  order  of  magnitude  thinner  than  the  core. 

Assuming  no  variation  of  displacement  through  the  thickness  the  exact  strain  displacement 
relationships  for  thick  circular  cylindrical  shells  are 


U;}  = 


d_ 

dx 


0 

0 


R  a 


R+z  dyj 

0  ; 


i 

R+z 


i  1  R  a 

R+z  1  R+z  dy 
I - 1 - - 

0  1 

_  _l 

<L  1 
dx 


dx 


~  .  0 


dx 


0 

0 


0 

zR  d 

R+z  dy 


zR  d 
R+z  dy 


R 

R+z 


0 

A. 

'  dx  _ 


w 

* 

cp 


'=[AJ  (5) 


_  R+z  dy 
where  =  u,  Uj  -  v,  etc. 

The  displacements  and  rotations  are  defined  in  Figure  47,  which  shows  the  coordinate 
system.  The  configuration  of  the  panel  is  shown  in  Figure  48.  Referring  to  this  figure. 
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FIGURE  47.  HONEYCOMB  SANDWICH  PANEL  COORDINATE  SYSTEM 
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FIGURE  48.  HONEYCOMB  SANDWICH  PANEL  CONFIGURATION 
AND  DIMENSIONS 
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the  panel  is  assumed  to  be  composed  of  four  individual  layers  of  material.  These  are: 
layer  1  -  the  core,  layer  2  -  the  inner  skin  (inner  refers  to  layer  closest  to  origin  of 
curvature),  layer  3  -  the  outer  skin,  and  layer  4  -  the  edge  attachment  or  "tabs."  In  the 
energy  analysis  a  presubscript,  r,  will  be  incorporated  to  indicate  the  layer.  This  notation 
will  be  necessary  for  the  strain  energy,  the  elastic  coefficient  matrix  and  the  strains. 


The  boundary  conditions  for  simply  supported  edges  are  assumed  to  be 

@  x  =  0  &  x  =  -t:  w  =  0,|i=0,|y=0 

@y  =  0&  y  =  b:  w  =  0,  |2  =  0,  |S  =  0 

The  constraint  on  w  is  a  pure  geometric  constraint;  however,  the  conditions  on  u,  v,  i|i, 
and  cp  are  necessary  to  satisfy  both  compatibility  relationships  and  geometric  relationships. 
The  du/^x  =  0  is  a  constraint  on  the  direct  strain  in  the  x  direction  and  satisfies  the 
moment  requirements  at  the  edge. 

Sinusoidal  mode  functions  are  chosen  as  the  assumed  mode  shapes.  They  are 


u  =  )  )  U  cosf3  x  siny  y 

Zj  Zj  mn  m  n' 

m=1  n=l 

v  =  )  )  V  sinP  x  cosy  y 

Zj  Zj  mn  m  '  n7 

m=l  n-1 

w  =  >  )  W  sinfJ  x  sin  y  y  > 

Zj  Zj  mn  m  n'  ' 

m=l  n-:l 

♦  =  I  I  ’.n  COsSmx  sinV 

m=l  n=l 

cp  =  )  )  $  si np  x  cosy  y 

^  L  L  mn  m  n7 

m-1  n=l 


(6) 


In  matrix  form, 


Cuk}=[Bj(q^ 


(7) 
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0 


0 


0 


B 


LW  0 

0  LXmY^J  0  0 

0  0  LXJfJ  0 

0  0  0  LXmYnJ 

0  0  0  0 


0 

0 

0 

WJ  J 


(7  cont'd) 


The  case  of  clamped  boundaries  will  also  be  considered,  provided  the  tapered  edge  is  not 
present.  The  clamped  edge  boundary  conditions  are 

w(0,y)  =  w(t,  y)  =  w(X/  0)  =  w(x,b)  =0 

wx(°,y>  =  wx(^/y)  =  wy(x,0)  =  wy(x,b)  =  o 

v{x,0)  =v(x,b)  =0 

cp(x,0)  =  <p(x,b)  =0 

♦  (0,y)  =t(t-,y)  =0 

u(0,y)  =  u{t,y)  =0 

The  clamped  edge  beam  functions  will  be  used.  These  functions  are  defined  in  Section 
II. A, 2,  equation  (6).  The  assumed  mode  functions  satisfy  the  notation  utilized  in  (7). 

Now,  substitution  of  (7)  into  (5)  gives 

te,j  =  [A]  [B]  (q^J  =  [ Dj  {q^n)  (8) 


and  substitution  of  (8)  into  (4)  results  In 


rUo  *  1  ^  C rD]T  [  C  ]  [  D  ]  (<£„)  -  $  Lq^J  [ ,E  Jlqkmn)  (9) 

where  the  presubscript  r  denotes  the  layer. 

The  strain  energy  of  the  vibrating  plate  is 

4 

u=  I  J,uodv  (10) 

r=l  vol. 

Integration  is  first  performed  with  respect  to  z  and  then  with  respect  to  x  and  y. 

The  integration  with  respect  to  z  is  represented  as  follows: 
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H/  / 


l2  b2^  h 


l\  bl 


-h, 


/  L^pqJ  [,0(ql}<l*+  /  Lq^J  C2E  3  Cqmn]d 

-h,  “bl”b2 


bl+b3 


+  /  L^pq-i  C3e3  Cqmn3<*a 


dydx 


\\  /'  f'Lq^^KqLj^^  /  / 


I  b  j  h^/2 


0  0  -H4/2 


0  b2-h4/2 
-t-  b2  h^/2 


b9  h i/2  „  , 

iJ  /  /  Lq^l  C4E  J^mnJd-dVdx  +  -5  /  /  /  Lqpqj[4E]iqmn}^ydx 


0  b,  -h4/2 

=  u1  +  u2+u3+u4+u5 


^2  ^1 


(11) 


The  first  three  terms,  symbolized  by  U.,  contain  the  strain  energy  of  the  core  and  face 
sheets,  whereas  the  last  Four  terms  represent  the  energy  of  the  tabs. 

The  types  of  integration  with  respect  to  z  which  appear  in  |E  are 


bl  ,  ,R  +  h.,  2h,  ,  h,  .  2h, 

/  ='"(Fr4)^('+^2)=-Hl 

-h. 


1 


h  .  /K  +  hi\  2hi 

-h 


R  +  hr 


J1  trl  =-2Rhi  -R2|n(trnr)^lT 


R  +  h  1  .  2  bl 


-h. 


/' 

-hr 


dz 


1  2hw,  ^_5ih 

HT'  R^vj""  r2  V  '  R2j  r2  H3 


(R+z)2  R+,,l 


(12) 
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*  (12  cont'd) 
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1 


-h 


1 


zdz 

(R  +  z)2 


R 

R+  h1 


R 

R  -  h] 


2 . 

z  dz  _  o.  K 

(R  +  z)2  ''R+hl 


R  -  hn 


The  required  integrations  in  £  are 


dz 

- 

J 

-hrh2 

R+z 

R  V 

-hi 

r 

dz 

J"2(\  - 

J 

_hrh2 

(R+z)2 

2  V 

R 

03) 


Integrations  appearing  in  3E  are 


J 


h,+h3 


dz  h3f,  2hl  +  h3  3h2  +  3hj h3  +  h2.  _ 
- 2T~  3r2  )  R 


r,+h3_*_JS(, 

J  R  +  z  p2  V  R  d2 


10 


04) 


The  indicated  integrations  over  the  surface  of  the  sandwich  portion  are  now  required.  The 
following  integrals  must  be  evaluated  to  write  the  E  matrices. 
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r  2 

X  X  dx  =  (1  -  6  ) 

J  pm  pm 


sin(Pc  ‘  em)t2  sin(Pp  +  + 

■  %  '  Pm>  2(VSm' 


sin(eD-BmK,,  sTn<ep  +  PmR, 

2(Sp  -  @m>  +  2(p„  +  PJ  J 


+  6 


pm 


p  '  m 

sinPmt2  COS  Pm't2  sinPmtl  COSpmtl 

2  "  2p  +  26 

m  m 


=  ,M 

1  pm 


. 


r  £ 

X’X1  dx  =  (1  -  6  ) 

J  p  m  pm 


sm(3D  -  3jt2  sin{3o  +  3mK-2 

%'^m>  VPJ 


sIn(3D-Pm^l  sin(gP4  PmHl 


+  6  3 

pm  m 


t-t2  -  -t,  si 

2  +- 


sin3m^2  co#ml2  singm^1  ^m*'! 


23 
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23 
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>05) 


=  0M 
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J  X" 


X  dx  =  -  3  iM 

pm  pi 
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p  l  pm 


X  X“dx  =  -  3^  i  M 
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'1 
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f  X"X"  dx  =  32  3^  ,M 
J  pm  p  m  I  pm 


Similarly, 
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I  Y  Y  dy  =  (1  -  6  n) 
J  q  n  7  q 


q  L 


sin(Yq  ‘  Yn)b2  s!n(Yq  +  Yn)b2^ 

*VYn>  ‘  2(YMq+Ynr 


sin(Yq-Yn)bT  ,  sin<Yg  +  Yn)bi  ,  . 

2(y  -  v  )  2(y  +  Y_) 

q  n  q  n  -* 

b2  -  bl  .  s'n^nb2  coiVnb2  ,  si"Ynbl  c°^nbl 


+  6 


qn  L- 


— 2 


*27 


2v 


=  ,N 

I  qn 


f  Y’Y'dy  =  (1  -  6  )y  Y 
J  q  n  7  '  qn  q  n 


sin(Ya  -  Vn)b2  sin(ya  +  V„)b2 

.  2^q  +  Yn> 


qn-  q 

sin(y  -  Yn)b,  sin(Y  +  Yn)b, 


v? 


%  +  Yn)  „ 


+  6  Y 

qn' 


•[ 


b2  -  bl  stnYnb2cosYnb2  sinY^  cosy^ 


2y. 


2y. 


(16) 


=  ,N 

2  qn 


f  Y"Y  dy  =  -  Y2  ,N 
J  q  n  7  q  1 


q  I  qn 


31 

b. 


f  Y  Y"dy  =  -  Y2  ,N 
J  q  n  7  nl  qn 


J1 

b. 


f  Y"Y"dy  =  Y2, 
J  q  n  '  q  I 


,N 

q  I  qn 


In  the  cqse  where  b]  -  0,  b2  =  b  and  t]  =0,  -t2  =  l,  the  above  integrals  for  simply- 
supported  edges  and  similar  integrals  for  clamped  edges  can  be  easily  written.  For  the 
case  of  the  constant  thickness  panel,  the  notation  of  (27)  and  (28)  in  Section  II  .A. 2 
should  be  used.  This  problem  has  been  individually  solved  and  a  report  is  being  prepared 
(Reference  16). 
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After  summing  the  first  three  integrals  in  (11)  and  performing  the  integrations  over  the 
surface,  the  integral  is  represented  as  follows: 


U1  =-5  /  /  Lq^pq-l  (iF+2F+3F)t<1mn3dydx"LqpqJ  (17) 


'1  D1 


where  K  is  a  generalized  stiffness  matrix. 


The  terms  in  the  matrix  are  given  in  Appendix  II. 

Now,  the  edge  numbers  will  be  treated  differently,  since  instead  of  being  thick,  as  in  the 
core,  they  can  be  considered  thin.  Therefore,  the  analysis  for  the  "tabs"  will  follow  that 
of  the  curved  plate  in  Section  H.A.2.  The  analysis  is  identical  up  to  the  point  where 
integration  over  the  surface  is  required. 


Therefore,  the  strain  energy  of  the  edges,  prior  to  integration  over  the  surface,  is  of  the 
form  of  equation  (16),  Section  II. A. 2.  integrating  over  the  surface  of  the  edges,  using 
terminology  of  equation  (16)  Section  II,  gives 


t  b. 


U 


l2  bl 


[ E ]  +  F] 


{q  Jdydx 
1  mnJ  ' 


The  matrix  of  coefficients  for  the  terms  in  5^  listed  in  Appendix  II. 

Integrations  of  the  mode  functions  over  the  surface  are  identical  to  the  terms  in  (15)  and 
(16)  except  for  a  change  in  limits.  Symbols  used  for  the  integrals  are  listed  below. 
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Limits  for  X  Integrations 
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Sp  <£  5Mpm 

Pp  A" 

Limits  for  Y  Integrations 


Integral 

bl  to  b2 

0  to  bj 

b2  to  b 

[  Y  Y  dy  = 

J  q  n  7 

,  N 

1  qn 

3  qn 

,N 

5  qn 

Jw 

2Nqn 

4  qn 

6Nqn 

/  YqYnd)'  ■ 

-  p2  ,N 
q  i  qn 

2 

-  9,  3Nqn 

-  p2  _N 
Kq5  qn 

f  Y  Y"dy  = 

J  q  n  7 

2 

-  r  n 

n  qn 

-  p2  -N 
n  3  qn 

2 

-  @n  5Nqn 

f  Y"Y"dy  = 

J  q  n  7 

p2  p2  ,N 
q  m  1  qn 

2  2 

P  p  nN 
q  n  3  qn 

p2  p2  SN 

Kq  n  5  qn 

r  i 

For  example,  the  integral  J  Y  Yndy  is  of  the  same  form  as  (16),  except  that  in  (16)  b j 

0 

is  replaced  by  bj  and  b^  is  replaced  by  0. 


The  integrals  tabulated  above  could  also  be  evaluated  for  clamped  edges;  however,  the 
size  of  the  task  was  prohibitive  for  this  investigation.  If  the  beam  functions  are  used. 
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each  of  the  above  integrals  would  contain  16  terms  to  be  evaluated  at  both  the  upper  and 
lower  limits. 


The  kinetic  energy  of  an  elastic  multilayered  shell  is 

T  =  I  /  {j  Pr(°r  +  \  +  *V)}dv 

r  vo  I . 


(19) 


where  r  indicates  the  layer  and  the  tir,  vr,  and  wr  are  total  velocities  of  a  point  (x,y,z) 
in  the  shell.  The  total  velocities,  defined  in  terms  of  mid-plane  velocity,  are 


Uj=U  +  Zl|j  ,  Vj  =  v  +  zcp  , 

u2  =  U  -  h|iir  ,  v2  =  v-h;cp  ,  w2 

U3  =  U+M'  v3  =  v  +  h1cp,  w3 


-  w 

=  w 

-  w 


u4  “  u 


v4  =  v 


w4=  W 


(20) 


Upon  introducing  the  actual  limits  of  integration,  (19)  becomes 

l2  b2  hl 

T  = -1  J  J  |p^  j  (u2+  2zii$  +  z^i|i2  +  v2  +  2zvqp  +  z2cp2  +  w2)dz  + 


l]  bl 


-hi 


-h. 


+  p2  /  (u2  -  2h1u$  +  h2^2  +  v2  -  2h1vcp  +  h2cp2  +  w2)dz  + 

-(h^h^ 


h1+h3 


+  p3  J  (u2  +  2h1u4i  +  h2ji2  +  v2  +  2h1vp  +  h2cp2  +  w2)dzjdydx  + 
bl 

l  b  h. 


(21) 


l  b^ 


+  -gp4/J  /  (u2+  v2+  w2)dzdydx  +  \  P4  J  /  /  (u2  +  v2  +  w2)dzdydx 

0  0  -h. 


^1  b2  b4 


0  b2  -h4 

t  b2  h4 


+  4  p 4  J  (  f  (u2  +  v2  +  w2)dzdydx  +  ^  P4  /  /  J  (u2  +  v2  +  w  )dzdyd? 


0  bj  -h4 


^2  bl  "b4 


Integrating  through  the  thickness,  introducing  the  assumed  mode  shapes,  and  integrating 
over  the  surface  yields  the  kinetic  energy  in  terms  of  generalized  coordinates. 


89 


(22) 


The  ;M's  ore  associated  with  the  different  surface  areas.  ]M  denotes  the  generalized  mass 
matrix  for  the  honeycomb  sandwich  center  sections  and  2^  f°  5^  are  ra I < zed  mass 

matrices  of  the  edges.  The  terms  in  the  M  matrices  appear  in  Appendix  II . 

It  is  of  interest  to  note  that  in  the  honeycomb  sandwich  analysis  the  generalized  mass  matrix  is  not 
diagonal  due  to  coupling  between  in-plane  velocity  and  shear  rotational  velocity. 

To  obtain  natural  frequencies  and  mode  shpaes,  harmonic  motion  is  assumed  and  Lagrange's 
equation  for  conservative  systems  is  utilized.  This  step  yields 


(23) 


This  is  the  form  of  the  eigenvalue  problem. 

In  applying  this  analysis  to  the  calculation  of  natural  frequencies  and  mode  shapes  of 
typical  structure,  it  is  necessary  to  define  the  elastic  coefficient  matrices.  For  a  honey 
comb  core,  the  elastic  coefficients  are 

]C11  =  1 C 1 2  =  1C21  =  1C22  =  0 


The  assumption  that  dw/^z  =  0  obviates  the  need  of  defining  ] C 1 3,  1C23,  1C32, 

and  1C33  since  they  do  not  appear  in  the  analysis.  In  actuality  Sw/5z  =  0  is  based  on  the 
assumption  that  the  core  is  inelastic  through  the  thickness  (the  C's  are  infinite). 


1 C .  j  —  1 G 

1  44  1  yz 

lC55=lGxz  ► 
1C66=  lGxy 


(24) 


The  facing  sheets  and  tabs  are  normally  constructed  from  isotropic  materials.  The  elasticity 
coefficients  are 


c , .  = 

E2 

r  - 

E3 

E4 

2^11 

1  -Uj 

3^11 

1  ~  u3 

4^11  ” 

- - 2 

1  “u4 

2C 1 2  =  2C21 

u2E2 

u3E3 

U4E4 

*  (25) 

'  ”  u2 

3^12“  3^21 

"  ,  2 

1  "  u3 

4C12  =  4C21 

E2 

r  - 

E3 

E4 

CM 

CN 

) 

CM 

.  2 

1  -u2 

3^22 

,  2 

1  -  U3 

4<-22 

CM 

D 
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2C66 


E2  E3  _  E4 

2(1  +  u2)  3C66  “  2(1  +  u3)  4C66  2(1  +  u4) 

(25  cont'd) 

When  this  theory  is  applied  to  practical  aircraft  structures,  the  following  simplifications 
can  be  made 


Hj  =  HjS*  1 

and  (h/R)Z  —  0  when  compared  to  1  in  the  same  expression.  When  used  separately,  this 
term  represents  the  contribution  of  bending  and  cannot  be  neglected. 


B.  Experimental 


1  .  Introduction 

This  section  covers  the  experimental  investigations  conducted  to  extend  the  analytical  and 
experimental  results  from  the  simple  panel  investigation  to  more  complex  flight  vehicle 
structure.  It  was  designed  to  determine  the  effects  of  dynamic  loads  and  heat  on  complex 
structure,  and  combined  the  environments  of 

a.  Low -frequency  vibration  of  the  type  encountered  during  taxi  or  high  speed 
buffet . 

b.  Aerodynamic  heating  typical  of  future  high  performance  aircraft. 

c.  Curvature  representation  of  typical  flight  vehicle  fuselage  structure. 

d.  Acoustic  excitation  from  broad-band  noise  source  typical  of  present  and  future 
power  plants . 

Titanium  honeycomb  sandwich  panels  were  selected  for  this  investigation  due  to  a  high 
strength-to-welght  ratio  and  good  physical  properties  in  the  temperature  region  involved. 

A  test  temperature  of  450°  F  was  selected  to  correspond  to  surface  temperatures  expected 
on  future  nigh  performance  vehicles. 

The  curved  honeycomb  sandwich  specimens  were  tested  under  a  combination  of  environment 
as  follows: 

a.  Broad-band  acoustical  noise. 

b.  Broad-band  acoustical  noise  and  heat. 

c.  Broad-band  acoustical  noise  and  low  frequency  vibration. 

d.  Broad-band  acoustical  noise,  heat,  and  low  frequency  vibration. 

The  objectives  of  this  phase  of  the  experimental  investigation  were: 

a.  To  determine  the  separate  and  combined  effects  of  low  frequency  vibratory 
loading,  heat,  curvature,  and  high  intensity  random  noise  on  the  fatigue  life 
of  titanium  honeycomb  sandwich  structure. 

b.  To  establish  design  criteria  for  sonic  fatigue  in  combined  environment. 
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c.  To  advance  the  state-of-the-art  for  sonic  fatigue  testing  and  life  prediction. 


2.  Test  Specimen  Design 

The  test  sDecimens  were  fabricated  of  titanium  faced  honeycomb  sandwich  to  meet  the 
requirements  of  the  thermal  environment  and  application  to  advanced  flight .vehicle  struc- 
ture  with  minimum  weight  and  fabrication  cost.  Some  fatigue  design  criterr 
Reference  17  were  used,  compensating  for  the  effects  of  curvature  and  elevated  tempera- 
ture  The  panels  were  designed  for  a  test  life  of  approximately  seven  hours  at  on  wefflN 
SPL  of  156  db,  five  hours  a?  an  overall  SPL  of  160  db,  and  one  hour  at  an  overall  SPL  of 
164  db,  giving  due  consideration  to  fatigue  data  scatter,  etc. 


Two  titanium  honeycomb  sandwich  configurations  were  designed.  The  details  are  shown  in 
the  following  table. 


TABLE  10 

TITANIUM  HONEYCOMB  SANDWICH  PANEL  DETAILS 


Item 

Pilot  Model 

Test  Specimen 

Facing  Sheets  (Ti,  6A1-4V) 

0.010  in. 

0.008  In. 

Doubler  (Ti,  6AI-4V) 

0.016  in. 

0.008  in.* 

Honeycomb  Core  (Polyester) 

3.5  Ib./fr 

3.5  Ib./fr 

Core  Thickness 

0.25  in. 

0.25  in. 

Adhesive  System 

HT  424 

HT  424 

Close-out  Member  (181  Fiberglass) 

2  Ply 

2  Ply 

Radius  of  Curvature 

84  in. 

84  in. 

Flat  Pattern  Size 

24  x  30  in. 

24  x  30  in. 

9 

Weight/Surface  Area 

- 

0.60  ib./fr 

*  0.020-inch  sheet  stock  chem-milled  to  0.016-inch  in  the  doubler  area 
and  0.008-inch  in  the  skin  area,  resulting  in  an  integral  skm-doubler 


arrangement , 


A  “pilot  model"  test  specimen  was  designed  and  built  to  evaluate  the  design  prior  to 
fabrication  of  the  fatigue  test  articles.  This  panel  was  tested  in  the  Lockheed  Hig 
Intensity  Sound  Facility  for  approximately  30  hours  of  broad-band  acoustical  noise  excita 
tion  at  160  db  overall  with  no  significant  failures.  Using  the  results  of  this  test  the  ina 
specimen  design  was  determined  to  be  that  which  is  shown  in  Table  10  and  Figure  49 
This  design  has  facing  sheet  thicknesses  which  are  considered  to  be  about  minimum  from 
the  fabrication  point-of-view . 
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3.  Test  Fixture  Design 


Test  Fixtures  were  designed  to  provide  peripheral  support  for  the  test  specimens 
representative  of  actual  flight  vehicle  installation.  The  design  required  the  capability  of 
providing  support  for  the  specimen  while  both  fixture  and  test  specimen  were  subjected  to 
low  frequency  vibration,  noise,  and  heat  environment.  Each  fixture  was  designed  to 
support  four  test  specimens  in  order  to  reduce  overall  test  time  and  cost,  A  limited  amount 
of  fixture  flexibility  was  required  to  introduce  shear  loads  into  the  specimen  when  the  low 
frequency  vibratory  load  was  applied,  A  detailed  stress  analysis  of  each  critical  component 
in  the  fixture  was  made  following  standard  aircraft  procedure  with  allowable  stresses 
obtained  from  Reference  18,  Each  part  was  designed  to  withstand  an  ultimate  alternating 
load  of  30,000  pounds  and  sound  pressure  levels  of  170  decibels.  A  typical  analysis  was 
comprised  of 

a.  Analysis  of  loads, 

b.  Determination  of  section  properties. 

c  ,  Tension  and/or  column  analysis , 

d  ,  Fastener  analysis . 

e.  Margin  of  safety  for  all  stresses. 

f.  Fatigue  life  at  all  stresses. 

The  fixtures  which  are  shown  in  Figures  50  through  52  were  capable  of  supporting  four  test 
specimens  simultaneously  in  a  manner  representative  of  actual  flight  vehicle  installation. 
The  upper  two  specimens  on  each  fixture  were  subjected  to  elevated  temperatures  produced 
by  seven  quartz  lamps  arranged  around  the  periphery  of  the  specimen  as  shown  in  Figure 
53,  Low  frequency  vibratory  loads  were  produced  by  an  electro-hydraulic  force  system. 

Test  specimens  were  installed  with  the  concave  side  toward  the  acoustical  noise  field,  or 
reversed  to  normal  flight  vehicle  installation,  to  protect  the  heat  lamps  from  the  high 
intensity  sound , 

4.  Test  Results 


a,  f reliminary  Tests .  Three  preliminary  tests  were  made  prior  to  the  fatigue  tests 
conducted  at  Wright-Patterson  Air  Force  Base,  Ohio,  These  tests  were  designed 
to  accomplish  three  things.  First,  a  "pilot  model11  test  was  used  to  evaluate  the 
fatigue  life  of  the  titanium  honeycomb  sandwich  design  prior  to  fabrication  of  the 
fatigue  test  articles.  Next,  tests  conducted  at  the  TF-33P-7  engine  test  stand  were 
designed  to  check  out  test  fixture  rigidity,  instrumentation,  and  quartz  lamp 
heating  system.  Last,  a  near-field  acoustical  noise  survey  was  made  to 
initially  set  the  position  of  each  test  fixture  in  its  desired  acoustical  noise 
environment. 

A  brief  description  of  each  preliminary  test  is  described  in  the  following  sections, 

•  "Pilot  Model1'  Test.  The  "pilot"  specimen  described  in  Table  10  was  tested  in 
the  Lockheed-Georgia  High  Intensity  Sound  Facility  prior  to  fabrication  of 
the  fatigue  test  specimens  used  for  this  investigation.  This  test  specimen  was 
instrumented  with  eight  uniaxial  strain  gages  strategically  located  to  define 
the  panel  response,  A  series  of  sine-sweeps  from  100  cps  to  700  cps  were 
made  to  determine  significant  frequencies,  damping,  etc.  A  plot  of  typical 
panel  strain  response  is  shown  in  Figure  54, 
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FIGURE  50.  COMPLEX  TEST  SPECIMEN  TEST  FIXTURE 


FIGURE  51.  FRONT  VIEW  -  TEST  FIXTURE 
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FIGURE  52.  FRONT  VIEW  -  TEST  FIXTURE  WITH  HYDRAULIC  ACTUATOR 


97 


■ALUMINUM  REFLECTOR  PLATE 


QUARTZ  HEATING  LAMPS 


REMOVABLE 

DUMMY 

INSPECTION 

PANEL 


FIGURE  53.  FRONT  VIEW  -  TEST  FIXTURE  WITH  TEST  PANELS  REMOVED 
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FIGURE  54.  TYPICAL  "PILOT  MODEL"  RESPONSE 


TABLE  1 1 

NEAR-FIELD  NOISE  SURVEY,  J57-P21  AT  MAX  A/B 


Microphone 

Position 

D  istance 
from  Engine 

Overall 
SPL  -  db 

1/3-0  BL 
db 

(160  cps*) 

1/3-OBL 

db 

(200  cps*) 

1/3-OBL 

db 

(250  cps*) 

1/3-OBL 

db 

(315  cps*) 

x/b 

Y/D 

i 

2 

0.85 

164.5 

152.5 

153 

153.5 

153 

2 

3 

1.03 

163 

146.5 

147 

148.5 

148.5 

3 

4 

1.20 

164 

147 

147 

148 

148 

4 

2 

1.56 

160 

147 

147 

148.5 

149 

5 

8 

1.91 

159 

147 

147 

148 

149 

6 

12 

2.61 

155.5 

145.5 

146 

144 

145 

7 

16 

3.32 

150 

138 

139 

138.5 

135 

8 

0 

2 

142 

121 

122.5 

124 

125.5 

9 

2 

2 

151 

131 .5 

132.5 

134.5 

135 

10 

4 

2 

161 

142.5 

142.5 

144 

145.5 

11 

4 

3 

147.5 

133 

133.5 

134 

137 

12 

6 

3 

155.5 

141 .5 

143.5 

143.5 

147 

13 

8 

3 

159 

145.5 

147 

146 

148.5 

14 

10 

3 

157 

145 

147 

146 

148 

15 

2 

4 

144 

128 

128.5 

126.5 

128.5 

16 

4 

4 

151 

134 

135 

134 

136.5 

17 

6 

4 

154.5 

139 

142.5 

142 

143 

18 

8 

4 

155 

141  .5 

145 

144 

144.5 

19 

10 

4 

156.5 

144 

147.5 

147.5 

146.5 

20 

12 

4 

156 

144 

148.5 

149 

147.5 

21 

5 

5 

512 

135 

138 

137.5 

138 

22 

10 

5 

154 

141 

146 

147 

143.5 

23 

15 

5 

151  .5 

139.5 

142 

144 

142 

24 

0 

10 

138 

116.5 

117 

120.5 

123.5 

25 

5 

10 

144.5 

120.5 

124 

128.5 

131 

26 

10 

10 

151 

128.5 

129 

135.5 

142 

27 

15 

10 

148 

133 

130.5 

135 

139 

‘Center  frequency  of  1/3  octave  band  filter. 
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Upon  completion  of  the  frequency  sweeps,  a  broad-band  noise  fatigue  test 
was  made  at  an  overall  sound  pressure  level  of  160  db  with  no  thermal  or  low 
frequency  vibration  input.  Thirty  hours  of  testing  at  160  db  overall  were 
accumulated  and  only  a  small  delamination  of  the  edge  member  was  the  result 
Data  obtained  during  this  test  were  used  to  establish  the  final  design  of  the 
test  specimen. 

•  Engine  Test  Stand  Checkout.  One  complete  test  fixture  comprising  an 
electro-hydraulic  force  system,  heating  lamps,  and  test  specimens  was 
checked  for  specimen  and  fixture  response  in  the  acoustical  noise  field  of  a 
TF-33P-7  jet  engine.  This  test  served  as  a  proving  ground  for  the  high 
temperature  strain  gages,  thermocouples,  and  adhesives  used  in  the  450°  F 
fatigue  test  thermal  environment. 

•  Near-Field  Acoustical  Noise  Survey.  Prior  to  positioning  the  fixtures  at  the 
test  site,  a  free-fleld  survey  of  the  acoustic  output  of  the  J57-P21  jet  engine 
was  made.  Sound  pressure  levels  were  sampled  at  62  positions,  27  of  which 
are  shown  in  Figure  55.  The  results  of  this  survey  are  tabulated  In  Table  11. 
Measurements  made  near  the  boundary  of  the  jet  efflux  were  used  to  position 
the  four  test  fixtures  in  the  desired  acoustical  noise  environment. 

b.  Fatigue  Tests  at  Wright  Field.  The  four  test  fixtures  were  placed  behind  the 
J57-P21  turbo- jet  engine  at  the  RTD  facility,  Wright  Field,  as  shown  in  Figures 
56  and  57.  Titanium  honeycomb  sandwich  specimens  1  through  16  were  tested 
during  the  first  part  of  the  test.  Heat  was  applied  to  each  of  the  upper  test  spec! 
mens.  Low  frequency  vibration  loads  were  applied  to  Fixture  2  and  3. 

The  tests  were  divided  into  segments  as  shown  in  Table  12.  Each  specimen  was 
thoroughly  inspected  after  each  segment. 


TABLE  12 

COMPLEX  SPECIMEN  TEST  RUNS 


Specimen  No. 

Segment  No. 

Run  Time 

Accumulated  Time 
Hr-Min . 

1-16 

1-6 

5 

00:30 

1-16 

7-9 

10 

01:00 

1-16 

10-17 

15 

03:00 

1-16 

18-23 

20 

05:00 

9-24 

24-51 

20 

14:20 

The  J57-P21  turbo-jet  engine  was  operated  at  afterburner  power  during  each  test 
segment.  It  developed  an  average  thrust  of  15,000  pounds  at  this  power  setting. 
The  quartz  lamp  heating  system  maintained  an  average  temperature  of  400°  F  on 
the  internal,  exposed  surface  of  each  heated  test  specimen.  The  electro- 
hydraulic  force  system  delivered  an  average  force  of  6000  lbs  RMS,  random 
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FIGURE  55.  NEAR  FIELD  NOISE  SURVEY-MICROPHONE  LOCATIONS 
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FIGURE  56.  OVERALL  VIEW  OF  TEST  ARRANGEMENT 
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FIGURE  57.  TEST  SET-UP 


amplitude  load  at  a  frequency  of  approximately  1  2  cps .  Strain  and  microphone 
data  were  recorded  on  magnetic  tape  during  each  run,  with  the  strain  from  the 
gages  mounted  at  the  center  of  the  test  specimen  as  shown  in  Figure  51  recorded 
continuously.  Microphone  data  and  strain  from  the  remaining  gages  were 
recorded  period tcolly  throughout  the  entire  test. 

Only  one  major  fatigue  failure  was  experienced  during  run  segments  1-24.  This 
failure  occurred  at  five  hours  of  test  time.  It  was  a  delamination  of  the  inner 
face  sheet  from  the  core  on  specimen  5,  Figure  58,  after  215  minutes  of  after¬ 
burner  time.  The  delmaination  was  deduced  to  be  the  result  of  a  void  in  the 
core-to-face  sheet  bond,  as  no  other  similar  failure  was  experienced  subsequently. 
Specimens  3  and  9  experienced  delamination  of  the  edge  member  around  the 
attachment  holes  as  shown  in  Figures  59  and  60.  Specimen  numbers  17  through 
24  were  then  installed  on  Fixture  1  and  2  in  place  of  specimens  1  through  8. 

The  low  frequency  vibratory  load  was  transferred  from  Fixture  2  to  Fixture  1  and 
testing  was  resumed  with  20  minute  time  segments. 

Testing  was  discontinued  after  accumulation  of  14  hours  and  20  minutes  of 
afterburner  time  with  only  nine  major  failures .  The  most  significant  fai  lure 
occurred  on  specimen  number  24,  Figure  61,  during  the  first  20  minute  run. 

It  was  caused  by  heat  and  a  void  in  the  bond  between  the  face  sheet  and  honey¬ 
comb  core . 

The  test  specimens  were  removed  from  the  fixtures  and  dye-checked  for  incipient 
cracks  at  the  completion  of  the  test  program.  Six  panels  had  skin  cracks  in  the 
inner  face  sheet  at  the  attachment  holes.  Figure  62  is  a  group  of  sketches  which 
show  the  location  of  each  crack  and  the  accumulated  time  at  the  end  of  the  test 

program . 

Microphone  and  strain  gage  data  recorded  during  the  test  were  analyzed  as 
fol  lows:' 

1  .  Microphone  data 

•  Overall  sound  pressure  level 

•  1/3-octave  band  sound  pressure  level 

2.  Strain  gage  data 

•  Root -mean-square  strain 

•  Narrow  band  10  cps  filter  bandwidth  (10-second  sample) 
m  Probability  density 

Typical  1/3-octave  band  sound  pressure  levels  are  presented  in  Figures  63  through 
65  for  specimens  1  through  5  and  10.  These  1/3-octave  band  analyses  of  the 
turbo-jet  noise  show  the  frequency  at  which  the  peak  sound  energy  occurs 
decreases  with  increasing  distance  from  the  exhaust  nozzle. 

Typical  narrow-band  strain  response  analyses  for  the  various  combinations  of 
environment  are  shown  in  Figures  66  through  69.  Tables  13  and  14  are  summaries 
of  the  measured  sound  pressure  and  strain  magnitudes  from  each  test  specimen,  as 
well  as  total  test  time. 
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SEPARATION  OF  EDGE 
MEMBER  FROM  INNER 
FACE  SHEET 


FIGURE  60,  PANEL  *9  (FRONT  VIEW)  TIME  TO  FAILURE  -  180  MINUTES 


SEPARATION  OF  OUTER 
FACE  SHEET  FROM  CORE 


OUTER  FACE  SHEET 


FIGURE  61. 


PANEL  #24  (FRONT  VIEW)  TIME  TO  FAILURE 


20  MINUTES 


107 


- - 

- - 1  I 

! ! 
'  I 
I 


f 


PANEL  #5 
TEST  TIME 
300  MIN. 


i 


I  I 
I  I 

I  I 

II 

I 


:tc- 


I 

I 

J.  I 

.1-SI 


ITT 


i 

i 

i 

i 

i 

i 

i 


PANEL  #  11 
TEST  TIME 
860  MIN. 


I 


I 


If 

"  "I- 

:i  r- 

:!  I 

'  ~1  I- 
f  II 

1  . 

■i  1 

it  1 

•i  i 

H  i 

PANEL  *  10 

•t 

TEST  TIME 

i  • 

•:i 

860  MIN. 

il: 

•i  1 
:i  i 

•1  u 

i  i  ■ 

1 1: 

1 1 f 

j 

•L 

.  j' 

± 


**  *  ■  X 

• n  - 

;i  r 

■  J  , 

i  i  ‘ 

*  1  l 

;t  i 
-1  i 

!! 

i  j 

SI  1 
*! 

i; 

PANEL  *14  1. 

TEST  TIME  . 

860  MIN. 

1  ' 

■i  i 

K  1 

i  : 

1 1  ■ 

i!L- 

1 1  ■ 
J  i; 

;  f - 

"  n 

■  i  i - 

■i  i 

j  i 

■  1 1 

1 1 
■  1 1 

i 

:i 1 

PANEL  *22 

! 

1 1  TEST  TIME 

I 

Tl  |  560  MIN. 

1 

■  f  I 

:i  i 

1 1 

•  i  i 

*i  i 

|i 

•1  L 

.j 1 

•}—  - - 

_ 1 

I 

I 


IM¬ 


PANEL  #23  |  t 

TEST  TIME  |  I 

560  MIN.  t  I 


!  I 

I  1 
l  L_ 

i — 


FIGURE  62  JEST  PANEL  SKIN  CRACKS 


108 


THIRD-OCTAVE  BAND  LEVEL  IN  DB  THIRD-OCTAVE  BAND  LEVEL  IN  DB 

(RE  0 . 0002  MICROBAR)  (RE  0.0002  MICROBAR) 


(a) 


(b) 


FIGURE  63.  MICROPHONE  DATA  -  1/3  OCTAVE  BAND  LEVELS 
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THIRD-OCTAVE  BAND  LEVEL  IN  DB  THIRD-OCTAVE  BAND  LEVEL  IN  DB 

(RE  0.0002  MICROBAR)  (RE  0.0002  MICROBAR) 


(o) 


(b) 


FIGURE  64.  MICROPHONE  DATA  -  1/3  OCTAVE  BAND  LEVELS 
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FIGURE  65.  MICROPHONE  DATA  -  1/3  OCTAVE  BAND  LEVELS 
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FIGURE  66.  COMPLEX  TEST  SPECIMEN  RESPONSE 
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FIGURE  67.  COMPLEX  TEST  SPECIMEN  RESPONSE 
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FREQUENCY  CPS 

FIGURE  68.  COMPLEX  TEST  SPECIMEN  RESPONSE 
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FIGURE  69.  COMPLEX  TEST  SPECIMEN  RESPONSE 


SPL  MEASURED  AT  PANEL  CENTERLINE 
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Note:  A  -  Sound 

B  -  Sound ,  Heat 

C  -  Sound,  Low-Frequency  Vibration 
D  “  Sound,  Heat,  Low-Frequency  Vibration 


TABLE  14 

STRAIN  MEASURED  AT  CENTERLINE  OF  CURVED  HONEYCOMB  SANDWICH  PANELS 


Frequency 

cps 

O  LO  o  o  o  o  o 

iO  O 

<NCS1CNCNCNCN(N 

270 

260 

1  1  1  o 

1  I  i  to 
CM 

1  o  o  o  o  o 

1  i  'O  K  to  u~y  sQ 
CM  CM  CM  CM  CM 

250 

Irt 

0) 

RMS  Strain 
H  in/in/cps 

<N(NO-tOOO 

'OfonMoO'Cs 

1  CM  o 

1  00  Os 

1  1  1  CM  | 
1  1  1  CO  1 

!  !  °  9  °  o 

e  i  in  —  ^  in  so 
i—  CM 

1  CM  1 

1  to  1 

u 

c 

ai 

□ 

cr 

p 

b> 

u_ 

HMi 

c 

o 

Frequency 

cps 

0  0  0  10  10  0*0 
*0  ''J‘  VO  *0  ''f  ''t  ■'O 

1  o  p 

1  O  K 

150 

1 

1  O  1 
t  O  l 

mZ 
.  C 

a> 

OO 

4- 

o 

c 

p 

RMS  Strain 

M.  in/in/cps 

O  O  O  MD  o  *0  o 
ococncoov^io 

CM 

1  COt  i 
1  —  CM  1 

i  1  i  OD  I 

!  1  1  CO  1 

1  COOOINO 
i  >o  rs.  io  ^  tj- 

100 

uo 

Frequency 

cps 

I  1  1  1  i —  ^ — 

i  <  1  1  r-  r- 

j  —  —  ! 

i  i  i  i 

iiii 

!  zz z z  | 

i  !  ! 

RMS  Strain 

H  In/in/cps 

1  1  1  1  o  o  o 

*  1  1  1  ON® 

O  ''t 

540 

500 

iiii 

(  i  i  i 

600 

750 

100 

150 

i  i  t 
i  i  i 

RMS  Strain 
|i  in/in/cps 

o  o  to  o  o  o  o 
io  o  o  o  o  o 

^tNCMM-CD  mo 

1600 

1250 

750 

1900 

2110 

725 

825 

1125 

1  o  • 

l  O  1 

o 

C 

0) 

E 

c 

e 

A 

A 

B 

B 

C 

C 

D 

D 

C 

C 

D 

D 

A 

A 

B 

flOUUQQ^  ^tflcQ 

£ 

LU 

c 

0) 

£ 

V 

Q> 

a 

o 

Z 

CMO^-tOONoOChO^  CNCO-M- 

117 


Note:  A  -  Sound 

B  -  Sound,  Heat 

C  -  Sound,  Low  Frequency  Vibration 
D  -  Sound,  Heat,  Low  Frequency  Vibration 


Figures  70  through  73  show  the  probability  distribution  of  strain  peaks  for  the 
various  combinations  of  environment.  It  can  be  seen  that  the  distribution  is 
essentially  Rayleigh  because  the  points  determined  from  experimental  data  lie 
near  the  Rayleigh  distribution  line  for  peak-to-rms  ratios  up  to  three. 


C.  Comparison  of  Theoretical  and  Experimental  Results 


As  a  means  of  confirming  the  theoreitcal  frequency  analysis  for  tapered  edge  sandwich 
panels,  one  of  the  titanium  panels  used  in  the  tests  described  in  Section  1H.B  was  tested  to 
determine  natural  frequencies. 


The  method  adopted  for  experimentally  determining  natural  frequencies  was  visual 
observation  of  Chladni  (nodal)  patterns  of  the  panel  when  excited  acoustically.  An  alumi¬ 
num  frame  with  overall  dimensions  of  24"  x  30"  and  a  cross-section  of  1"  x  4"  was  used  to 
clamp  the  edges  of  the  panel .  This  frame-panel  combination  made  up  the  front  wall  (grill- 
cloth  side)  of  a  dual  speaker  enclosure.  The  surface  of  the  panel  was  four  inches  from  the 
two  15"  Altec  605A  loudspeakers.  The  speakers  were  electrically  connected  so  that  the 
phasing  could  be  reversed.  This  allowed  determination  of  some  of  the  even  modes. 


The  edge  frame  was  rigid  enough  to  stabilize  the  panel;  however,  it  did  not  produce 
clamped  edge  boundaries. 


The  experimental  edge  condition  was  probably  closer  to  simply-supported  than  clamped  as 
verified  by  the  test  comparison. 

Panel  dimensions,  physical  constants,  and  elastic  constants  were 
=  0.75",  l2  =  22.00",  t  =  22.75" 


b1  =  0.75",  b2  =  28.00",  b  =  28.75", 

h|  -  0.100",  h2  =  0.008",  h3  =  0.008",  h4  =  0.020" 

P1  =  5.25  x  10  6  #-sec2/m4,  P2  =  P3  =  P4  =  4.15  x  10~4  #-sec2/in4 
=  1 1 ,000  psi,  jGyz  =  21 ,500  psi 

E2  =  E3  =  E4  =  1  *6  x  *°7  Psi'  v2  =  v3  =  v4  “  °-322 

The  facing  sheets  were  6AI-4V  titanium  and  the  core  was  polyester  honeycomb.  The  edge 
member  was  a  laminate  of  6AI-4V  titanium  and  fiberglass.  For  calculation  purposes,  it 
was  considered  to  be  titanium  alone  because  the  fiberglass  does  not  contribute  significantly 
to  the  total  stiffness. 

The  following  table  shows  calculated  and  measured  natural  frequencies. 
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PROBABILITY  OF  PEAKS  EXCEEDING 


FIGURE  70-  PROBABILITY  DISTRIBUTION  OF  STRAIN  PEAKS 
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PROBABILITY  OF  PEAKS  EXCEEDING 


FIGURE  71.  PROBABILITY  DISTRIBUTION  OF  STRAIN  PEAKS 
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PROBABILITY  OF  PEAKS  EXCEEDING 


FIGURE  72.  PROBABILITY  DISTRIBUTION  OF  STRAIN  PEAKS 
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PROBABILITY  OF  PEAKS  EXCEEDING 


FIGURE  73.  PROBABILITY  DISTRIBUTION  OF  STRAIN  PEAKS 
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TABLE  15 

NATURAL  FREQUENCIES  OF  TITANIUM  SANDWICH  PANEL 


n 

m 

i 

2 

3 

i 

244.0 

204 

399.0 

347 

593.1 

617 

Calculated 

Experimental 

2 

219.0 

219 

412.6 

642.6 

627 

Calculated 

Experimental 

3 

364.5 

369 

521  .4 
484 

750.9 

Calculated 

Experimental 

Although  the  theory-experimental  agreement  is  good  for  a  majority  of  the  modes,  it  is 
believed  that  differences  in  the  higher  modes  are  due  to  inaccuracies  in  the  transverse 
shear  moduli  for  the  honeycomb  core.  In  the  next  section,  the  effect  of  transverse  shear 
modulus  on  natural  frequency  is  shown. 

Also,  as  stated  earlier,  the  actual  boundary  conditions  of  the  test  panel  were  not  the  same 
as  the  boundary  conditions  of  the  theory.  This  is  the  most  probable  cause  of  the  differences 
between  theoretical  and  experimental  values  in  the  lower  modes. 


D  .  Discussion  of  Results 


1 .  Theoretical  Results 

A  number  of  parameter  variation  studies  were  made  to  assess  the  importance  and  the  effects 
of  individual  parameters  on  the  natural  frequencies.  To  keep  the  studies  simple,  only  one 
parameter  was  varied  for  each  study.  All  variations  were  made  around  the  basic  pane! 
design  used  for  comparison  of  theory  and  experiment-  Some  parameters  such  as  core  thick¬ 
ness  and  radius  are  extremely  influential  in  the  determination  of  frequency,  whereas  others 
such  as  tab  width  have  only  secondary  effects.  The  effect  of  these  parameters  an  natural 
frequency  will  be  discussed  in  the  following  section. 

a.  Core  Transverse  Shear  Modulus  of  Elasticity.  In  general,  experimental  values  of 
transverse  shear  modulus  have  not  been  very  reliable.  It  is  probably  the  most 
difficult  parameter  to  accurately  measure  of  all  the  parameters  associated  with  a 
honeycomb  sandwich  design.  For  this  reason,  the  effects  of  this  parameter  are 
considered  to  be  very  important. 

Two  parametric  studies  were  made  in  assessing  the  effects  of  shear  modulus.  In 
one  case,  the  two  shear  moduli,  Gxz  and  GyZ,  were  set  equal  and  varied  over 
the  range  of  100  to  1,000,000  psu  The  variation  of  modal  frequency  with  shear 
modulus  is  shown  in  Figure  74. 

One  very  important  conclusion  gained  from  this  investigation  is  that  core  shear 
modulus  has  only  a  minor  effect  on  the  natural  frequency  of  the  first  mode.  It 
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Transverse  Shear  Modulus  of  Elasticity  in  Psi 
FIGURE  74.  EFFECT  OF  VARIATION  OF  TRANSVERSE  SHEAR  MODULUS 


serves  only  as  a  means  to  separate  the  two  facing  sheets.  This  implies,  that  for 
practical  purposes,  theories  which  define  first  resonance,  neglecting  transverse 
core  shear,  will  give  acceptable  results.  This  is  substantiated  for  flat  panels  by 
Sweers  (Reference  19). 

The  effect  of  transverse  shear  modulus  on  high  frequency  modes  is  appreciable. 

In  this  case,  low  values  of  transverse  shear  modulus  tend  to  cause  a  decrease  in 
resonant  frequencies.  The  frequencies  are  considerably  lower  than  those  calcu¬ 
lated  by  a  theory  which  neglects  the  transverse  shear  modulus.  In  Reference  19, 
the  following  simple  equation  for  calculating  natural  frequency  is  given: 


where 


(2h1  +  h2)zh2E 

D  = - — 

2(1  -/) 


Higher  mode  frequencies  for  values  of  transverse  shear  modulus  above  100,000  psi 
could  be  calculated  using  the  simple  frequency  analysis  shown  above.  In 
instances  where  the  core  has  a  weak  transverse  shear  modulus  the  theory  needs 
refining . 

Reduction  of  the  basic  curved  design  used  in  these  parameter  studies  to  a  flat 
panel  gives  a  calculated  1,1  mode  of  92.6  cps  and  721  cps  for  the  3,3  mode 
frequency.  The  simple  theory  gives  106  cps  for  the  1, 1  mode  and  954  for  the 
3,3  mode.  Although  this  is  a  relatively  thin  panel  with  high  shear  modulus,  the 
effect  can  still  be  observed.  The  ratio  of  simple  theory  to  complex  theory  for  the 
1st  mode  is  1.15  while  It  increased  to  I  .32  for  the  3,3  mode.  If,  instead  of  the 
high  strength  fiberglass  core,  a  core  with  shear  modulus  of  1,000  psi  had  been 
used,  the  above  comparison  would  show  an  error  factor  of  approximately  two. 

Another  interesting  effect  is  observed  by  holding  Gxz  constant  at  11,000  psi  (the 
basic  design  value)  and  varying  GyZ  from  2000  to  500,000  psi.  Figure  75.  The 
results  are  similar  to  Figure  74.  Tne  investigation  emphasizes  the  importance  of 
core  weakness  and  shows  that  a  core  which  is  weak  along  one  transverse  axis  will 
have  lower  natural  frequencies  in  the  higher  modes  than  an  isotropic  core. 

b.  Effect  of  Curvature.  Curvature  is  the  most  significant  parameter  from  the  stand- 
point  of  its  effect  on  frequency.  In  Figure  76,  the  radius  is  varied  over  a  broad 
range  about  the  basic  panel  design  radius  of  84  inches,  at  a  radius  of  9  inches 
and  an  arc  length  of  27.25  inches,  the  basic  panel  design  will  subtend  an  angle 
of  180  degrees.  The  radius-to-thickness  ratio  in  this  case  is  greater  than  40; 
therefore,  the  approximation  of  plane  sections  remaining  plane  after  deformation 
is  valid . 


At  a  radius  of  1,000",  the  panel  is  flat  for  all  practical  purposes.  As  expected, 
the  lowest  frequency  mode  is  the  1,1  or  fundamental .  At  the  basic  panel  design 
radius  of  84  inches,  the  1,2  mode  has  the  lowest  frequency  with  the  fundamental 
modal  frequency  being  slightly  higher.  This  "cross-over"  effect  was  noted  by 
Arnold  and  Warburton  in  a  theoretical  study  of  the  natural  frequencies  of  thin 
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FIGURE  75.  EFFECT  OF  ONE-DIMENSIONAL  VARIATION 
OF  TRANSVERSE  SHEAR  MODULUS 
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FIGURE  76.  EFFECT  OF  CURVATURE 
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cylinder  (Reference  20),  The  experimental  results  do  not  demonstrate  the 
"cross-over"  predicted  by  theory.  The  1,1  mode  was  measured  at  204  cps  and 
the  1,2  mode  was  219  cps  so  the  trend  to  cross-over  is  evident  in  the  results. 

Examination  of  the  curves  in  Figure  76  reveals  that  modal  frequencies  which  have 
node  lines  parallel  to  the  straight  edges  of  the  panel  (1,2  and  1,3)  are  insensitive 
to  a  change  in  radius  as  compared  to  modal  frequencies  with  node  lines  parallel 
to  the  curved  edges  (2,  !  and  3,  l), 

c,  Effect  of  Core  Thickness.  The  simple  theory  for  flat  panels  included  in  the 

discussion  of  shear  modulus  is  of  the  following  form,  provided  core  thickness  is 
the  only  variable: 


f  =  C  (2h|  +  h9) 

n  n  1  l 

This  relationship  Is  plotted  in  Figure  77  for  the  first  mode  frequency.  The 
constant  C  was  determined  using  the  base  panel  first  modal  frequency.  In  com¬ 
paring  this  simple  relationship  to  the  thickness  effects  calculated  with  the  theory 
of  Section  III  ,A,  it  is  seen  that  the  above  equation  does  not  provide  a  first  order 
approximation  of  the  actual  results  for  the  lower  order  modes.  The  slope  of  the 
higher  modes  does  approach  that  of  the  linear  approximation  for  an  intermediate 
range  of  core  thickness. 

Actually,  the  first  mode  frequency  variation  shown  in  Figure  77  should  probably 
be  attributed  to  curvature  rather  than  thickness  for  thickness  less  than  0.3", 

This  phenomenonhas  also  been  noted  in  frequency  calculations  of  thinshel  Is  where 
it  was  observed  that  the  "cross-over"  was  more  predominate  for  thinner  shells. 

It  is  interesting  to  note  that  higher  modes,  which  were  not  affected  by  curvature, 
for  Instance  the  1,3  mode,  have  approximated  the  same  slope  as  that  of  the 
linear  relationship.  The  slope  at  higher  thickness  is  attributed  to  transverse  shear 
effects , 

d.  Effect  of  Skin  Thickness.  The  express  purpose  of  the  core  in  a  honeycomb 

sandwich  is  to  stabilize  the  facing  sheets  In  creating  a  rigid,  light-weight  struc¬ 
ture.  In  designing  the  honeycomb  sandwich  panel  it  is  desirable  to  create  the 
greatest  stiffness-to-weight  ratio  possible  tn  order  to  decrease  the  bending 
stresses.  Therefore,  an  Item  for  assessing  a  honeycomb  panel  design  is  its 
stiffness-to-weight  ratio.  In  the  simple  flat  honeycomb  panel  theory,  the  natural 
frequency  is  proportional  to  7D/p  *  For  very  thin  faces,  where  the  core  controls 
panel  weight,  the  frequency  varies  approximately  as  the  square  root  of  face 
thickness  .  In  the  mid-range,  the  frequency  is  nearly  independent  of  Face  thick¬ 
ness  and  when  skin  thickness  approaches  the  order  of  magnitude  of  core  thickness, 
the  frequency  approaches  a  linear  variation. 

The  effects  described  above  are  somewhat  different  for  the  curved  honeycomb 
sandwich  panel.  Examination  of  the  curves  in  Figure  78  shows  that  frequency  at 
first  increases  with  skin  thickness,  then  peaks  and  subsequently  decreases  for 
increases  in  facing  thickness.  It  is  the  peak  in  frequency  that  is  of  significance 
in  structural  design.  The  peak  in  frequency  represents  the  most  efficient  structure 
from  the  standpoint  of  rigidity. 
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FIGURE  77.  EFFECT  OF  CORE  THICKNESS 
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FIGURE  78.  EFFECT  OF  SKIN  THICKNESS 


Another  interesting  aspect  of  the  effect  of  face  thickness  is  concerned  with  the 
ratio  of  inner  skin  to  outer  skin  thickness.  In  Figure  79,  it  is  seen  that  if  both 
skins  are  of  approximately  the  same  thickness  the  natural  frequencies  are  highest 
indicating  the  highest  ’’dynamic  stiffness.”  This  is  also  indicative  of  efficiency 
in  design. 


2.  Experimental  Results 

a.  Test  Specimen.  The  test  specimen  for  the  complex  specimen  investigation  was 
designed  using  an  engineering  approach.  State-of-the-art  design  criteria  were 
taken  from  Reference  17  to  initially  set  the  dimensions.  The  design  nomographs 
from  Reference  17  were  modified  as  described  in  Appendix  I  to  account  for 
curvature  effects  and  the  effect  of  elevated  temperature  on  fatigue  life.  Extrapo¬ 
lation  of  the  nomograph  of  Reference  17  was  necessary  to  obtain  a  design. 

A  "pilot  model"  of  the  design  was  made  and  tested  in  the  High  Intensity  Sound 
System  Facility  to  experimentally  determine  the  fatigue  life  of  the  design.  It 
was  exposed  to  30  hours  of  broad -band  acoustical  noise  at  an  overall  sound  pres¬ 
sure  level  of  160  decibels  without  a  significant  failure.  The  "pilot  model"  design 
was  too  strong  so  it  was  weakened  using  the  strain  data,  etc.,  from  the  test  to 
scale  it  down.  As  a  result,  the  facing  sheets  were  reduced  from  0.01"  to  0.008" 
and  the  edge  doubler  was  also  reduced  from  0.01"  to  0.008"  .  The  total  metal 
thickness  at  the  edge  was  0.016"  .  From  the  production  point  of  view,  most  of 
the  design  details  were  minimal .  The  test  specimen  had  an  average  surface 
weight  of  0.6  pound  per  square  foot,  which  Is  less  than  most  structural  components 
designed  to  withstand  environmental  conditions  comparable  to  those  of  this  test. 

b.  Test  Fixture.  The  test  fixtures  were  designed  to  provide  peripheral  support  for  the 
test  specimens  representative  of  flight  vehicle  installation.  In  addition,  fixture 
flexibility  was  necessary  to  transmit  bending  and  shear  loads  produced  by  the 
electro-hydraulic  shaker  to  the  honeycomb  sandwich  panels.  The  test  fixtures 
were  designed  to  hold  four  specimens  at  one  time  so  that  each  environmental 
condition  could  be  evaluated  simultaneously  to  minimize  the  test  time  and  cost. 

During  the  fatigue  tests  at  Wright  Field,  Ohio,  a  few  of  the  fixture  test  specimen 
support  structures  experienced  fatigue  cracks.  These  support  members  were 
replaced  with  a  like  part  and  the  tests  were  continued.  Stop  drilling  was  also 
used  but  only  as  a  temporary  measure.  The  cracks  are  believed  to  have  been 
caused  by  a  horizontal  component  of  the  low-frequency  vibratory  force  resulting 
from  lateral  vibration  of  the  test  fixtures.  This  lateral  vibration  was  produced  by 
Jet  efflux  buffet  and  was  not  accounted  for  in  the  stress  analysis. 

c.  Acoustical  Noise  Environment.  Figures  63  through  65  show  the  1 /3-octave  band 
sound  pressure  levels  in  the  range  of  significant  test  specimen  response  was  about 
140  decibels.  The  maximum  acoustical  energy,  in  most  cases,  appeared  at 
frequencies  above  315  cps.  This  is  well  above  the  frequencies  demonstrating  the 
greatest  test  specimen  strain  response.  Downstream  from  the  engine  nozzle. 
Fixtures  3  and  4  positions,  the  acoustical  energy  peaked  at  lower  frequencies, 

250  cps.  Jet  efflux  buffet  of  the  fixture  was  an  undesirable  feature  of  these 
locations,  however. 

A  comparison  of  the  sound  pressure  levels  from  the  near-field  noise  survey 
contained  in  Table  1 1  with  those  made  at  the  surface  of  the  test  specimens,  Table 


131 


1000 


sd^  'Xouanbaj-j  juouosay 


132 


FIGURE  79.  EFFECT  OF  VARYING  RATIO  OF  INNER  SKIN  THICKNESS  TO  OUTER  SKIN  THICKNESS 


13,  reveals  that  the  near  field  survey  1/3-octave  band  levels  were  higher.  The 
difference  is  attributed  to  three  things: 

•  Condenser  microphones  were  used  for  the  near-field  noise  survey  and 
piezoelectric  microphones  were  used  far  test  monitoring  purposes. 

•  Ambient  temperatures  were  greater  for  the  fatigue  tests  than  for  the  near-field 
noise  survey  so  the  acoustical  noise  output  of  the  J-57-P21  turbojet  engine 
was  greater  during  the  near-field  noise  survey. 

•  The  near-field  noise  survey  was  made  in  the  absence  of  reflective  surfaces, 
except  the  concrete  ground  plane,  whereas  the  test  fixtures  provided  reflective 
surfaces  during  the  fatigue  tests. 

d.  Test  Specimen  Strain  Response.  A  study  of  Figures  66  through  69  shows  four 
distinct  strain  response  patterns.  Figure  66,  the  specimen  response  to  acoustical 
noise,  shows  the  greatest  response  occurred  at  160  cps  or  the  (1, 1)  mode  with 
some  activity  in  the  (2,1)  mode,  250  cps,  and  the  (1,3)  mode,  300  cps.  When 
heat  was  applied  to  the  test  specimen  the  strain  response  was  diminished  in  the 
(1,1)  mode,  (2,1)  mode,  and  (1,3)  mode,  as  shown  in  Figure  67.  Also,  the 
frequency  of  each  mode  was  lowered  significantly  at  the  400°  F  temperature  as 
was  the  case  for  the  elevated  temperature  sonic  fatigue  tests  of  the  simple  titanium 
panels  described  in  Section  II.B.  Figure  68  shows  the  response  of  the  low  frequency 
random  vibration  and  acoustical  excitation  was  the  response  of  the  test  specimen 

to  the  low  frequency  random  excitation  superimposed  on  the  response  due  to 
acoustical  excitation.  The  strain  was  quite  large  at  12  cps  and  the  (1,1)  mode 
was  about  the  same  as  that  for  acoustical  excitation  alone.  Effect  of  all  three 
environments  is  shown  in  Figure  69.  The  strain  response  had  the  same  general 
appearance  as  that  shown  on  Figure  67  except  the  strain  at  12  cps  is  considerably 
higher  because  of  the  low  frequency  random  loading. 

The  higher  modes  such  as  (2,1)  and  (1,3)  had  magnitudes  almost  equal  to  and  in 
one  case  greater  than  the  magnitude  of  the  (1 , 1)  mode.  The  higher  modes  were 
influenced  by  the  unequal  distribution  of  acoustical  energy  in  that  their  excitation 
was  3  to  5  decibels  more  than  the  (1,1)  mode  excitation. 

e.  Probability  Distribution  of  Strain  Peaks.  Probability  distribution  of  strain  peaks 
was  determined  from  strain  data  recorded  during  the  fatigue  tests.  Figures  70 
through  73  are  the  probability  distribution  plots  for  each  of  the  environmental 
conditions.  It  can  be  deduced  from  these  probability  distribution  plots  that  the 
strain  peaks  for  all  environments  had  a  distribution  which  differed  only  slightly 
from  a  Rayleigh.  The  strain  peaks  for  the  higher  modes  did  not  differ  appreciably 
from  the  fundamental  mode,  and  the  strain  peak  distribution  at  12  cps  was  very 
nearly  Rayleigh. 

f .  Fatigue  Test  ■  There  are  four  possible  modes  of  fatigue  failure  for  honeycomb 
sandwich .  Cracks  may  form  at  the  fastener  holes  around  the  edge  of  the  panel, 
cracks  may  form  in  the  facing  sheets,  the  honeycomb  core  may  experience  shear 
fatigue,  and  the  facing  sheet  may  separate  from  the  core  (bond  failure). 

The  test  specimens  experienced  two  of  the  four  types:  core-to-facing  sheet  bond 
failure  and  cracks  at  fastener  holes.  The  most  common  type  of  failure  experienced 
in  flight  vehicle  service,  facing  sheet  cracks,  was  not  produced  even  though  the 
specimens  were  subjected  to  extremely  severe  test  conditions.  Only  nine  of  the 
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twenty -four  test  specimens  were  damaged.  If  fatigue  damage  is  used  as  a 
criterion  for  measuring  the  severity  of  each  environmental  condition,  the 
following  categories  of  severity  can  be  devised.  First,  the  most  severe  environ¬ 
ment,  according  to  the  number  of  fatigue  failures  and  time-to-failure,  is  heat 
and  acoustical  excitation.  Three  of  the  nine  failures  took  place  under  these 
conditions.  Next,  the  environment  of  low  frequency  random  vibration  and 
acoustical  excitation  should  be  placed  in  second  position  even  though  it  also 
produced  three  damaged  test  specimens.  (The  average  time-to-failure  was  con¬ 
siderably  more  for  this  test  condition  than  the  average  time-to-failure  for  the 
environment  of  heat  and  acoustical  excitation.)  Acoustical  excitation  alone 
produced  two  damaged  test  specimens  so  this  environmental  condition  should  be 
placed  in  third  position.  Finally,  one  test  specimen  experienced  fatigue  cracks 
when  exposed  to  the  three  environmental  factors  simultaneously. 

Before  the  tests  were  conducted,  it  was  believed  that  the  environment  comprised 
of  heat,  acoustical  excitation,  and  low  frequency  random  vibration  would  be  the 
most  severe.  No  explanation  can  be  offered,  but  there  is  some  evidence  that  the 
combined  environment  is  less  severe  than  the  other  environmental  conditions.  For  exam¬ 
ple,  test  specimens  7  and  8,  Figure  57,  were  mounted  in  Fixture  ^3  during  the 
first  part  of  the  test  program  and  were  exposed  to  the  combined  environment  for 
over  300  minutes  without  failure.  These  specimens  were  removed  and  test  speci¬ 
mens  23  and  24  were  mounted  in  their  place.  The  electro-hydraulic  force  system 
was  also  removed  and  then  attached  to  Fixture  .  Consequently,  test  specimens 
23  and  24  were  subjected  to  heat  and  acoustical  excitation  only.  After  20 
minutes  of  testing,  one  of  the  facing  sheets  of  test  specimen  23  separated  from 
the  honeycomb  core  producing  a  core  bond  failure.  When  the  final  inspection 
was  made,  test  specimen  24  had  cracks  around  the  fastener  holes. 

It  was  not  possible  to  correlate  time-to-failure  with  the  environmental 
conditions,  so  sonic  fatigue  tolerance  curves  are  not  included.  The  honeycomb 
sandwich  design,  however,  looks  very  promising  for  use  In  design  of  future  air¬ 
craft.  Due  to  its  light  weight,  approximately  0.6  pound  per  square  foot,  and 
apparent  high  resistance  to  sonic  fatigue,  it  would  be  desirable  for  high  perfor¬ 
mance  flight  vehicle  structure. 
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IV 


CONCLUSIONS 


1  .  The  curved  panel  frequency  analysis  theory  derived  herein  yields  accurate  results, 
especially  for  the  higher  modes. 

2.  Curvature  greatly  increases  the  sonic  fatigue  resistance  of  simple  panels. 

3.  The  sonic  fatigue  design  nomograph  for  skin-rib  construction  which  accounts  for 
structural  curvature  effects  will  increase  sonic  fatigue  design  capabilities  and  aid  in 
saving  flight  vehicle  weight. 

4  Sonic  fatigue  tests  mode  on  8AI-IM0-IV  titanium  alloy  and  17-7PH  stainless  steel 
flat  panels  demonstrated  that  exposure  to  heat  and  acoustical  excitation  simul¬ 
taneously  was  many  times  more  severe  than  exposure  to  heat  and  then  acoustical 
excitation  alternately. 

5.  Stresses  resulting  from  acoustical  excitation  and  thermal  buckling  were  the  primary 
cause  of  the  fatigue  of  the  fatigue  cracks  developed  during  the  elevated  temperature 
tests  of  the  simple  panels,  not  a  degradation  of  fatigue  properties  due  to  elevated 
temperatures.  (See  Reference  21  .) 

6.  The  frequency  theory  for  cuived,  tape  red -edge,  honeycomb  sandwich  panels  derived 
herein  is  an  accurate  first  step  toward  improving  honeycomb  sandwich  sonic  fatigue 
design  techniques. 

7  The  complex  specimen  test  program  was  an  undertaking  that  was  much  too  ambitious. 
It  should  have  been  conducted  in  several  phases,  investigating  one  environment  at  a 
time  and  then  all  possible  combinations. 

8.  A  jet  engine  is  a  poor  acoustical  noise  generator  for  conducting  sonic  fatigue 
research  investigations  because  there  is  no  way  of  controlling  the  spectrum. 

9  An  analytical  and  experimental  program  should  be  devised  to  increase  the  reliability 
and  range  of  application  of  the  existing  honeycomb  sandwich  design  nomographs. 

10.  The  complex  specimen  tests  showed  that  lightweight  honeycomb  sandwich  can  be 

designed  to  withstand  environmental  conditions  expected  to  be  associated  with  high 
performance  flight  vehicles  of  the  future. 
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APPENDIX  I 


DESIGN  NOMOGRAM  MODIFICATIONS 

The  effect  of  curvature  on  stress  ratio,  expressed  by  equation  (53),  Section  II  .A. 3,  can  be 
used  in  conjunction  with  existing  design  charts  for  flat,  rib-skin  structure.  The  theoretical 
values  of  the  constants  C  and  D  in  (53)  do  not,  however,  yield  sufficiently  accurate  results 
to  be  applied  to  a  design  analysis.  If  the  results  of  tests  described  in  Section  II  .B  are  used 
to  calculate  the  constants,  then  a  design  chart  can  be  produced  empirically. 

The  constants  C  and  D  have  the  following  average  values  determined  experimentally: 

C  =  0.006 

D  =  0.412 

Using  the  above  values  (53)  is  presented  as  a  nomogram  in  Figure  1-1  .  Using  a 
conventional  design  nomogram  for  flat  panels  such  as  that  shown  in  Figure  t-2  (from 
Reference  17)  an  estimation  of  panel  dimensions  can  be  obtained.  Using  these  dimensions, 
a  stress  reduction  ratio  due  to  curvature  is  calculated  from  Figure  1-1 .  The  procedure  is 
repeated  using  the  material  fatigue  curves  for  reduced  stress.  The  calculation  is  iterative 
and  2  to  3  iterations  should  be  sufficient  for  satisfactory  convergence. 
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FIGURE  1-2  .  DESIGN  CHART  SKIN  AND  RIB  CONSTRUCTION 


APPENDIX  II 


FINAL  EQUATIONS 


The  elements  of  the  K  matrices  are  listed  below. 
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